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ABSTRACT! Automatic Digital Encoding System II (ADES II) is a system 
for the automatic translation of mathematical formulas into programs 
of coded instructions for an electronic digital computer. ADES 
represents a new approach to the problem of automatic programming in 
that it is applicable to almost all modern computers, and it is not 
limited to a special class of mathematical problems. The system 
consists of a formulation language, an Encoder, and a digital computer. 
The formulation language closely resembles ordinary mathematical 
language, and is based on the theory of recursive functions. The 
Encoder is an automatic device which receives a mathematical formulation 
as input and produces the necessary computer instructions as output. 
The logical design of an Encoder is given. It consists of a computer 
with certain routines loaded into its storage. Certain desirable 
features of the computer in such a system are suggested. 
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This report contains the results of research directed toward expediting 
the process of preparing, or programming, mathematical problems for 
computation on high-speed electronic digital computers. The advent of 
the automatic computer which performs computations at electronic speed, 
has revolutionized the art of computation, but it has raised problems 
of its own. For many mathematical problems, the speed of computation 
is so great and the process of programming so involved, that the time 
spent on the computer in useful computation is often Just a fraction 
of the total time elapsed between the formulation of the problem and 
the obtainment of the numerical results. 

This report presents a system designed to correct this situation. The 
system relieves the mathematician of a major portion of the programming 
task. In effect, it replaces the human by an automatic device which 
carries out the programming duties. The design for such a device is 
given. An experimental model is in process of construction. 
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PREFACE 



At the September 1955 Annual Meeting of the Association for 
Computing Machinery, the author presented a paper which described 
Automatic Digital Encoding System I, (ADES I), a system of automatic 
translation of mathematical formulas into computer programs. The 
present paper describes an extension of that system, ADES II, which 
significantly enlarges the class of problems which can be automatically 
programmed. It is conjectured that ADES II is sufficiently general to - 
cope with most of the two-dimensional mathematical problems that are 
submitted to existing digital computers. 

ADES I excluded problems which involved double recursion. This 
is now provided for in ADES II, which is also capable of programming a 
special kind of triple recursion. However, ADES II harbors certain 
obvious limitations which can and should be removed. Further, there is 
room for much improvement in procedure, for example, in input-output 
procedure. Finally, there remains the deeper problem of efficiency,both 
in time and storage utilization. The minimalization or efficiency problem 
has thus far been treated as secondary in urgency. These considerations 
motivate plans for another system, ADES III, now in the research stages. 

The ADES approach to automatic programming is believed to be entirely 
new. Mathematically, it has its foundations in the bedrock of the theory 
of recursive functions. The proposal to apply this theory to automatic 
programming was first made by C. C. Elgot, a former colleague of the 
author *s. While at the Naval Ordnance Laboratory, Elgot did some research 
on a language for automatic programming. Some of his ideas were adapted to 
ADES, and we wish to acknowledge this fundamental contribution. 

It is our belief that the theory of recursive functions, as propounded 
in references QlJ , [2\ , and [5J say, and the closely related mathematical 
logic provide the natural, possibly- the only, framework for constructing 
a general automatic coding system for digital computation. However, it 
appears that these subjects are not yet within the purview of the majority 
of those mathematicians presently concerned with preparing problems for an 
electronic digital computer. Since this report is directed primarily at 
them, ADES H is presented here as a practical system for the translation 
of mathematical problems Into computer programs, rather than as an 
application derived from the theory of metamathematics. The mathematical 
logician will find no theorems to justify the methods adopted, and the 
sporadic references to recursive function theory are intentionally in an 
intuitive vein. The style is consciously expository and heuristic rather 
than telegraphic. However we do assume that the reader has a good mathe- 
matical background. ADES 18 designed primarily for the mathematician. 
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This report is in two parts, under separate covers. The first part 
presents definitions and rules for formulating problems in a stylized 
mathematical language. The second part describes the logical design of 
the 'Encoder*, a machine which will receive the mathematical formulation 
as input and yield a complete computer program as output. The logical 
mechanism of the Encoder becomes rather intricate in places. Therefore, 
minute details are often omitted, since it is felt that the mathematician 
versed in machine techniques will be able to supply the missing details, 
in one form or another, once he has grasped the overall pattern. In fact, 
he can modify the design in many ways, for there usually is more than one 
means to an end. Many features in the present design were suggested by 
the author's colleague, Leroy Krider, Mr, Krider and the author are working 
on the construction of an actual Encoder for ADES II, which will utilize a 
650 Magnetic Drum Calculator, A Joint report describing this Encoder and 
experiences with its operation will be published in the near future. 

Many thanks are due to Larry Schmid for his critical reading of the 
manuscript and helpful suggestions. We also thank Mrs, M, Zawatzky for 
her assistance in preparing the manuscript for publication. 
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AUTOMATIC DIGITAL ENCODING SYSTEM II (ADES II) 

Introduction , In the last few years, we have witnessed the advent of 
automatic computing machines which perform digital computation at 
electronic speed. In fact, the speed of computation is so great that 
for many mathematical problems, the time spent on the machine is hut a 
small fraction of the total time elapsed between the formulation of the 
problem and the obtainment of the numerical results. 

Much of the time and effort required to solve any moderately complex 
mathematical problem by means of an automatic digital computer is consumed 
in the process commonly referred to as 'programming'. 

Programming can be described as that process which starts with the 
problem as a set of statements in mathematical language (i.e. equations), 
coupled with a set of metamathematical statements in ordinary language 
(e.g. English), and transforms these two sets into a single set of instruc- 
tions coded in a language intelligible to the machine. This set of coded 
instructions, known as the 'program*, is a list of the successive computer 
operations required to obtain the numerical results. 

To elaborate, programming can be considered to take place in two stages. 
The first stage involves the writing of explicit mathematical formulas which 
define the quantities to be computed, indicating the order in which these 
quantities are to be computed, $..e. the flow chart), and finally, specifying 
the numerical values to be assigned to the independent variables in these 
formulas. The writing of the formulas may be a matter of transcription, or 
it may entail some numerical analysis and algebra, for example, converting 
implicit equations into explicit ones, setting up a numerical integration, 
or an interpolation, and so on. 

In the second stage of programming, the mathematical formulas obtained 
in the first stage, together with the metamathematical statements which 
specify the flow chart and the procedure for substituting numbers for the 
variables, are translated into the language of the machine. This translation 
is essentially a clerical process in which three main tasks are performed. 
(1) The computation is broken down into an ordered sequence of machine 
operations. These include arithmetic operations and logical operations, 
each of which is translated into a numerical code. (2) The numbers to be 
substituted for the independent variables, (i.e. the input data) are 
assigned storage locations, and their addresses are then combined with the 
appropriate operation codes from step (1) to form machine instructions. 
(5) Storage locations are allocated to the computed results, and machine 
instructions for the output of these results are given. 
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These three tasks of encoding are recognized to be quasi-mechanical, 
and, to a great extent, they can be assigned to a machine, preferably to 
the computer itself. Toward this end, various schemes such as compilers, 
library routines, interpretive routines, and pseudo-code systems have been 
devised. As yet, these schemes have fallen short of the immediate 
objective, which is to accept a mathematical formulation of a problem as 
obtained In the first stage of programming and produce a complete program 
for machine computation. 

In the present paper, a system is presented, which, for the most part, 
realizes this objective for a wide class of problems. This system will be 
called *Automatic Digital Encoding System, II* , (ADES,II). The Roman 
numeral, II, indicates both the incompleteness of the system and its 
extensibility. The incompleteness is of two kinds. First, certain 
components of the system are still in an unrefined though workable state, 
resulting in a loss of efficiency. Second, a class of mathematical problems 
cannot be encoded by system II, namely, those problems which necessarily 
involve arbitrary triple recursions or, more generally, recursions on k 
indices, where k > 5. However, the system does provide for a special 
£Lnd of triple recursion which arises frequently in practice. Furthermore, 
there is no inherent reason for excluding the general triple recursion. 
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1. The Components 

As conceived in this paper, an automatic encoding system consists 
of three primary components! a formulation language, an Encoder (i.e. 
a machine which does the encoding), and an electronic digital computer. 
The three components influence each other to such an extent that it is 
difficult to construct any one of them independently of the others. 
Seemingly, this is a disadvantage, since in most instances we are given 
the computer and required to design the other components around it. It 
would appear that each computer requires a system peculiar to itself. 
Actually, the computer does not play a critical role. 

Fortunately, there are certain general principles, drawn from the 
theory of metamathematics, which are applicable to all digital computation. 
Thus, although ABES II has been set up for computers with certain character- 
istics, the system can be modified very easily to fit most digital computers. 
In particular, the formulation language and the logical design of the 
Encoder can be used universally with relatively minor changes, while the 
internal machinery of the Encoder can readily be adapted to different 
computers . 

2. The Computer 

To fix our ideas, we assume that the computer in ADES H is a 
single-address, sequenced-program, floating-point calculator which can 
perform the usual arithmetic and logical operations. It can also carry 
out a new operation which facilitates the metamathematical processes of 
the computer, i.e. those processes which compute and modify storage 
addresses. This operation is called •Modify*. 

The coded instruction, *Modify n», causes the next instruction in 
sequence to be executed as if the number in storage n had been added to 
it. For example, if the number 230 is in storage n, then the pair of 
instructions, Modify n. Multiply 1500, is equivalent in effect to the 
single instruction, Multiply 1550, i.e. multiply by the number in 1550. 
Since n is arbitrary, this makes it possible to use any storage in the 
computer as a •B■*ox , , or 'index register 1 . Furthermore, the instruction, 
•Modify n», converts the floating-point number in storage n into an 
integer suitably positioned for address modification. This eliminates the 
need for two kinds of arithmetic in the computer, and makes the mathematical 
formulation of a problem conceptually simpler. It also simplifies the 
machinery of the Encoder. 
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The above features are not absolutely essential in system II. 
However, they are convenient. In any case, they can be achieved by means 
of interpretive routines if the actual computer hardware is lacking. 

5. Requirements of the Language 

The formulation language is the critical component in any automatic 
coding system. In devising such a language, one is beset by conflicting 
demands. The language should resemble 'ordinary mathematical notation* 
as much as possible. It should not result in inefficient programs. It 
should not be so flexible that the Encoder becomes excessively complicated. 
Its alphabet and syntax should be rich enough to accomodate any effectively 
calculable function. 

The first three properties depend largely on personal Judgment. In 
the desi gning of system II, there was a tendency to choose a la ng ua g e 
which simplified the Encoder at the expense of efficiency in the computer 
program. The resulting inelegancies can be removed, and this will be done 
In system HI. Most Important is the fourth requirement, and to satisfy it. 
recourse was had to the theory of recursive functions (see t,l3 > X 2 J » t'l '• 
This theory is pertinent, since the arithmetic and logical operations 
performed by a digital computer are, in effect, operations on integers, and 
what's more, they are 'primitive recursive' functions. 

In the domain of positive integers, one distinguishes two kinds of 
functions, primitive recursive and general recursive. Intuitively speaking, 
a primitive recursive function is one which is obtained by composition of 
a finite sequence of arithmetic functions and recursions on one index. A 
general recursive function is obtained by adducing multiple recursions on 
k indices, where k can be as large as we please. Theoretically, multiple 
recursions, as such, can be eliminated if a new operator, ji, called the 
minimization operator, is included in the l an g ua g e, [lj . Thus, 

4 (x^ ..., x n )= ^y( f(Xj., ..., x n , y) * o} ) 

means that <J(x., ..., x^) is the mini mm value of y such that f^, ...,x n ,y) = 0. 

The minimization operator appears to be of limited practical value in 
eliminating multiple recursions. Hence, the language for system II provides 
for explicit double recursions, and a special type of triple recursion, as 
explained later. 

In the next few sections, the language will be described send-formally 
hy presenting its alphabet and syntax in the form of definitions and rules, 
and illustrating them with numerous examples. 
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k. The Alphabet and Elementary Expressions 

The alphabet is constructed from the nine generic symbols , 
q,r,a,b<,c,f,e,p and d, the positive integers and zero. The usage of 
the generic symbols is indicated by their names. 

q... .Independent index 

r.„.. .dependent index 

a.... independent variable (or data symbol) 

b. dependent variable (or defined variable) 

c... ............. .free variable 

f. function symbol 

e...... equal sign, ( = ) 

p ••••punctuation symbol 

d.. output symbol 

The letters of the alphabet consist of the generic symbols with 
numerical subscripts attached. Thus, the alphabet consists of the 
independent indices q^q,,^, ..., the dependent indices r ,r ,r~, ..., 

the independent variables a.a,,a 2 ,..., and so on for the other symbols. 

In what follows, we shall abbreviate, using •variable 1 to denote 
either one of the independent variables a ,a,,ap, ..., or one of the 

dependent variables b ,b.,,b , ..., unless otherwise stated. Likewise, 

'index 1 will mean either an independent or dependent index. 

Def. fr-.l . An indexed variable of degree one is an expression of the 
form, xi, where x denotes a variable and i denotes an index* 

e.g. I ^V b^, a^r^, b^ . 

An indexed variable of degree two is an expression of the form, xij, 
where x denotes a variable, and i and j denote indices. 

e.g. : ajO^qg, b^g^, b^r^, a^rg . 

An expression consisting simply of a variable is said to be an 
indexed variable of degree aero , 

e.g. : a^ b 2 , a v * 

A numerical constant is a floating-point number, written in some 
standard form consisting~bf an exponent, a modulus, and algebraic sign, 
(Note I The precise form cannot be specified, since it depends on the 
computer. However, in this discussion, the exact form is irrelevant, and 
we shall write numerical constants in the form that is easiest to present.) 
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Definition 4.1 can be extended to indexed variables of 
arbitrary degree. However, in ADES, II, the maximum degree of a variable 
is two. 

In the formulation of a problem, the independent variables are used 
to denote the various kinds of input dataj e.g. a might represent 
pressure, a^ temperature, a time, a. the elements°of a given matrix, etc. 

If a, is to assume n numerical values as the computation proceeds, the n 
values being supplied as a row of data, then in the formulation this data 
must be represented as an indexed independent variable of degree onej 
e.g. a-o^, where q. ■=■ 0,1,2, ..., n-1. (Thus, in that formulation, a. 
must always occur with one index.) For data supplied in the form of a 
two-dimensional matrix, we must use an independent variable of degree two J 
e.g. a^g-q,, where q^ and q, assume successive positive integral values 

between specified lower and upper bounds. A data constant must be denoted 
either by an indexed variable of degree zero, say a,-, or by its numerical 
value. 

A computed quantity (i.e. an intermediate or final result) which is 
defined by a designated formula is represented by a dependent variable. 
If several values of a dependent variable, b_, are to be computed, they can 
be distinguished, if necessary*, by using one or two indices J e.g. b 1 q 1 ,b_q_g <) . 

It should be apparent from these remarks that the independent indices 
in a problem are to assume successive integral values as the computation 
proceeds . The lower and upper bounds of each independent index used must 
be specified explicitly. This specification is called 'quantification' of 
the index, and is an important part of the formulation. A quantification 
consists of a phrase 'for all integral values of q,, such that L,^ q,.fS u, ». 
In this report, we shall use the quantification symbol, * \/ 'to represent 
•for all*, and a quantification will be written briefly in a form to suggest 
• V <3h * ^h 6 li ^ "i ' • Here, L, and u, denote independent variables or 

indices which specify the bounds of q.. Rigorous rules for quantification 
will be given in section 8. 

It is often necessary to have indices which do not take on successive 
values, or it may be convenient to have an index which is defined in terms 
of a quantified index. In such cases, dependent indices are used. For 
example, if a problem calls for an index which takes on only even integral 
values, we define r^ 2 q,, where q.ss 0,1,2,... ,n, and write a-r instead 
of a. (2q 1 ) . Other uses or the dependent indices will appear later. 

5. No tation 

In ADES, the formulas which define the dependent variables and 
dependent indices must be written in what is known as parenthesis—free , 
or prefixed-operator notation, m this concise notation, the formula 

*The necessity depends on other factors to be explained later. 
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W"*" a 2^ ** wri " tten as + a i a o* and a i' a 2. ** written *» • a i a o« ^ other 
words, each operator (i.e. function) symbol is written to tEe left of its 

operands, and all parentheses are eliminated since they are redundant. 

Note that juxtaposition, &••& , does not connote multiplication. The 

centered dot is used for tne^multiplication operator. 

For a unary operator, such as cosine, the prefixed-operator notation 
coincides with ordinary mathematical notation, provided that all parentheses 
are eliminated, e.g. cos x, A unary operator is said to be of degree 1, a 
binary operator (e.g. ■♦■,•) is of degree 2, and in general, a function of 
n operands is of degree n. In the actual ADES II language, each function 
is denoted by the symbol f with identifying subscript. For example, f 
denotes +> f, denotes — , f ^ denotes • , and so on. However, to simplify 
the exposition, in this report we shall continue to use the conventional 
mathematical designations +-, • , cos etc., with a few exceptions. One of 
these is the identity function, which is usually not written explicitly in 
conventional mathematics. In ABBS, it will be denoted by f , and must be 
written wherever called for by the syntactical rules. Thus, as we shall 
see, one must write b = ^a, instead of b. » a. . Note that for expository 
reasons again, we shall write ' -= ' for an ordinary equal sign instead of 
•e *, as required in ADES. Likewise, in this paper, all common punctuation 
in°?ormulas will be written in conventional form instead of using the 
symbol P with a subscript. 

Now, to further illustrate the use of prefixed-operator notation, we 
list several examples written in both conventional and prefixed-operator 
notation. 

Conventional P refixed-Operator 

0*1+ a 2 ) a^.. 

/fel<*L>-«fe< r l* P 2fl- 

(%+ qjK*! - r 2 ) 

sin (d^ a^)) 
2.^ + 3.) 

6. Well-formed formulas and terms . 

An expression in parenthesis— free notation consisting of a juxtaposition 
of operators and operands is called a 'string 1 DO . Not all strings are 
meaningful J e.g. +ap, sin a-b., +■ JLaa^ are meaningless. A meaningful 

string is called a 'welJU-formed formula*, (w.f.f.) . In essence, a string 
is a well-formed formula if for each operator in the string there is the 
correct number of operands. A rigorous definition is as follows. 



• +a l a 2 a 3 




<*-V^p Cfc— **Q 




V- w* 


: r l r 2 


••HL^- 


■ r l r 2 


sin + 


b 1 a 1 q 1 


•2.4 


• 8^3. 
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If x is an operand, then x is a vel^forrned formal*. If y L ,y 2 „.„y n 

are welk-formed formulas and f. is a function of degree n, then 
f b y_y 2 . # .y is a well-formed formula . 

In AXES, we use 'terms' rather than well— formed formulas. 

Def. 6.1 , An index term is a well-formed formula in which the operands 
are independent indices, q^q^qg, ..., dependent indices, r^r^r^..., 

indexed independent variables, or numerical constants, and in which there is 
at least one function symbol « 

Hotel The w.f ,f„, a., is not an index term, nor is the w.f,f., q. • 
Instead, one must write f.a- and f-q. , where f. is the identity function, 

Def, 6,2. An index equation (or r-equatlon ), is an expression of the 
form 

where J denotes one of the dependent indices r > r i> r 2 > •••» and 4 denotes 
an index term. (Note the comma after (J.) 

We illustrate this definition by examples written both in conventional 
and in ADES language. 



Conventional ADES 

•f* o « 

9 



r l= \ iq 2 h r l- f l 3 l q 2' 

r o-*L-V r o Sm Vl.> 

r 5 = fc^) + r 2 ] / q^ r 5 =, /+ a^r^, 

r 6 = q 1 (q 1 +a 1 ) / /2. , r g - / - ^ + ^2., 



Def. 6.5 . A b— term is a well— formed formula in which the operands are 
indexed variables or numerical constants, and in which there is at least one 
function symbol. 

Def. 6,4. A b-equation is an expression of the form, 

y — * > , 

where y denotes one of the dependent variables h ,b.,b , ..., and t denotes 
a b— term. 



8 
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e.g. 

Conventional ABES 



b i = [ V*l> + a a3 a 3 \ - - + W*V 

b 17" b 5 (r l' r 2 ) - a l ( V q 2 ) b 17~ ^l'sVlV 

It may be necessary to have an independent Index, q, say, as an 
operand in a b-equation. This is done by introducing an independent 
variable, a.., and defining a^q.) » q . The identification of a i with 

q is done in the ^Computer Table*, to be explained in section 13. 

J 

Def. 6.5. An integer-valued independent variable is an independent 
variable which assumes the values of a designated independent index. 
Integer-valued variables are permissible operands in b-equations. 

The b-equations and the r-equations constitute the major part of the 
formulation of a problem, since they define what is to be computed. 
However, there is need for another kind of equation. Suppose a formula 
occurs several times in the same problem. It would sometimes be convenient 
to be able to write this formula once with 'free* or unspecified variables. 
Then each time the formula is used, it would be necessary to indicate only 
which data are to be substituted for the free variables. For this purpose, 
we adduce the *free terms * and Auxiliary equations* defined as follows. 

Def. 6.6 . A free term is a well^-formed formula in which the operands 
are free variables or constants, and in which there is at least one function 
symbol. 

Def, 6.7 . An auxiliary equation is an expression of the form, 

where g denotes a function and () a free term. The degree of g is equal to 
the number, n, of free variables in <f. If x^ ..., x n denote permissible 

operands for an index term, then gx,, ..., x is the index term obtained 

when x. is substituted for the free variable c. in <| ; and similarly for 
b-terms. The function, g, on the left side or"an auxiliary equation is 
called an «auxiliary* function. In this report, we shall use subscripts 
of 50 or greater to designate auxiliary functions. 
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For example, suppose the index equations, 
r ± - <L X (% f 1), r 2 ^ qgCg^-J- 1.), ^ = q^ + 1.) , 

occur in the same problem. If we define an auxiliary function, f say, 
by the auxiliary equation, - 30 

f 50 = * c l + c 2 1# ' 
then we can write the index equations as 

r l~ f 50^1^-' 
r 2 = f 50«2V 

Note the following important point, which we emphasize as a Rule. 

Rule 6.1 . Each equation in an ABES formulation must be punctuated 
by a comma at the right end of the equation. 

The functions used in an ADES formulation are of two types^ the 
auxiliary functions defined above, and •library* functions. To explain 
what is meant by •library function 1 , we must refer to the computer again. 

7. Functions and the Computer . 

Associated with a digital computer is a library of subroutines for 
computing various mathematical functions, m ADES, it is assumed that 
the computer library contains at least the floating-point arithmetic 
operations, + ,-, • ,</. In some machines, these subroutines are part 
of the hardware. In others, interpretive subroutines are used. In either 
case, we regard them as part of the library of functions, that is, as 
mathematical functions which can be called for without special formulation. 

Each library function is designated by an f with a characterising 
subscript, (in ADES II, we have reserved subscripts to ^9 for library 
functions. As mentioned earlier, f, is the identity function, f_ is ■+* , 
f is -, f^ is • , t c is/ , and so on.) In the translation process, the 

Encoder recognizes the subscript and inserts the pertinent calling sequence 
for the library subroutine into the program. At the same time, it compiles 
the subroutine itself, that is, it arranges the transfer of the subroutine 
from external storage (e.g. tape) to a suitable location in high speed 
storage, (see description of Encoder.) 

The main point here is that in formulating a problem, a library 
function can be called for simply by writing an f with the correct subscript. 
Ho special equations are needed to define the function.*-*- On the other 



*1. See Appendix I for exceptions. 
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hand, each auxiliary function must be defined by an auxiliary equation. 
The assignment of a subscript to an auxiliary function is arbitrary, 
provided that it does not conflict with the subscript of a library function 
or another auxiliary function in the same problem. 

We assume that the ADES II computer library contains, besides the 
four arithmetic operations, the following elementary functions I absolute 
value, square root, direct and inverse circular and hyperbolic functions, 
exponential, logarithm, power function. 

In the terminology of recursive function theory, the library functions 
are the given functions of our formal system. The formulas for b's, r»s 
and auxiliary f *s are derived from the library functions by applying the 
syntactical rules. These rules are such that any derived function is 
primitive recursive relative to the set of library functions. Conversely, 
the library and rules of syntax should be such that any primitive recursive 
function can be derived. This implies that the syntax should include 
definition by recursion. Recursion will be introduced in section 11. In 
ADES II, however, recursive definitions are restricted to dependent variables, 
that is, a dependent index cannot be defined by recursion. Thus, not every 
primitive recursive function can be formulated. However, in practice, only 
a narrow class of functions is excluded thereby. 

8. Quantification, Phasing J 

The b-equations constitute the main part of an ADES formulation. They 
define the quantities to be computed, and closely resemble the conventional 
equations which a mathematician writes to describe a computational problem. 
But, in automatic digital computation, mathematical formulas are not 
sufficient. One must also specify the order in which quantities are to be 
computed, that is, a flow chart must be drawn up. This is frequently more 
difficult than writing the formulas. 

In ADES, a large part of the flow chart is determined implicitly and 
in a natural way by the very structure of the b-equations. However, certain 
explicit directions must be written J e.g. for the quantification of 
independent indices, for branch equations, and for recursions. This section 
presents some of the rules for writing quantifications. 

Def. 8.1 . a quantifier is an expression of the form * V Liu*, where 
the symbol » V •* is the universal quantifier 'for all', i denotes an 
independent index, L denotes the lower bound of i and may be either an 
indexed independent variable or an index, and u denotes the upper bound 
of i and may be an indexed independent variable or an index. ~* 

Thus, * \/Liu* is to be read as 'for all integral values of i such 
that L^i$u.« 

* In the actual ADES alphabet, Y is denoted by a library function symbol. 
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Rule 8 ' 1 * A quantifier is written to the left of a b-equation. 
Unless specifically forbidden by other rules, any b-equation may be 
preceded by one or more quantifiers. (The b-equation and the dependent 
variable which it defines are then also said to be quantified.) If more 
than one quantifier is written to the left of a b-eguation, the order of 
quantification is considered to take place from left-to-right, i.e. if 
we have 

v h.h\ v Va u 2 y - 4* 

this Implies that for all i^, L ± * i^S t^, y is to be computed for all 

i 2 ,L 2 4 i 2 ^ Up* The , I aantificati0,^ ' V L l i x u x , apples to the 
quantification • V LgigUg 1 • Hence, L g and u g can depend on 1^, whereas 
L. and u. must be independent of ig. 

It is sometimes necessary to quantify several b-equations by the same 
quantifier, or it may be necessary to specify the order in which certain 
independent quantifications are to take place, m many problems, it is 
necessary that certain b*s be computed before others. These situations 
axe provided for by the 'phase equation' . 

Def. 8.2 . A phase equation is a b-equation of the form, 

y=>Vi y 2 ••• V » 

where the y*s denote dependent variables, n is an arbitrary positive integer, 
and f is a special library function . 

The function, f , is to be read as, 'compute the following quantities 1 . 
Strictly speaking, f° is not a function in the mathematical sense. However, 
in ADES, it is to be°regarded as a function of arbitrary degree which 
operates on those b's between itself and the first comma to its right. 
(Therefore,the explicit and correct placement of this comma is essential, 
Just as it is in all b-equations.) 

The dependent variable denoted by y on the left side of a phase equation 
does not represent a computed mathematical quantity; but, this introduction 
of a dependent variable allows us to compound phase equations, for y can 
be used as an operand in another phase equation, or even in an ordinary 
b-equation. This makes it possible to specify the order of computation by 
the simple expedient of writing the dependent variables in the desired 
order within a term. This is amplified in Appendix H, and In examples. 

Rule 8.2 . To quantify several dependent variables *>$!»•••> b j n > 
with the same quantifier, a phase equation is written as follows: 

YLiu b jQ = f Q Vj2 — V # 
(This can be read as «for all i, L 4 i £ u, compute To^t •••> b jn # *' 
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Def. 8.3 . A phase equation in which the left member is the special 
dependent variable, b o , is called a master phase equation; i.e. 

b z f b.,b._ ... b. . • 
o o 31 32 Jn' 

Rule 8.3 . There is precisely one master phase equation in each problem 
formulation. The dependent variable, b , is used only on the left side 
of the master phase equation and in no 8ther equation. 

The master phase equation is the first b-equation to be programmed. 
The left-to-right order of the b.,, ..., b. in the master phase equation 
determines the order in which thi various pnases of the computation will 
proceed. We shall postpone a precise detailed description of the ordering 
procedure to Appendix II, but many of the main ideas are suggested in the 
following examples, and a rule of thumb is given in example 8.4- below. 



e.g. 8.1 . It is required to compute a 1 (q :L )+ agCq-J, and a^q^ - a^q.^ 

for all q,,0 4 q.4 25 . A formulation, neglecting input-output symbols, 
is as follows* 



Voq^ 



b o * f o V>1 ' 
bj* 4-a L q 1 a 2 q 1 , 

b 2 s -a^a^ , 



The master phase equation specifies that b. and b_ are to be computed for 

all q 1 ,0 ^^4 25. The order in which the b*s on the right side of a 

b-equation are written, reading from left to right, determines the order 
in which they are to be computed. In all problems, the computation starts 
with the master phase equation. In example 8.1, the quantification informa- 
tion will be obtained first. Then the programs for b and b_ will be set 
up in that order. Note that, in this case, the order in which the b-equations 
is written is irrelevant. This is true for all b-equations except recursion 
equations and vector equations (see sections 10,11). 

e * g * 8 - 2 ' "W* :» * ****** *° CamSttlte *** *L<*L> + V^ aad 
*$(<l 1 ,<3 2 )+ ajjq^qg), for all q^O 4 q-^ 9, and for all q^O 4 <3g4 8 » 
A formulation, neglecting input-output symbols, is as follows! 

b l = 4" V^lW 
^ b 2 - + «j^^ 
VOq^ VOqgS b o = f^bg, . 
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Starting with the master phase equation, the program for the quantification 
of g^ will be set up first. The quantifier for cu will be set up next, and 
then the programs for b and b in that order. Again, the b-equations can 
be written in any order. 

e.g. 8.3 . It is required to compute 3}^)+ agC^) for all ?i' ^ <&£ 5* 
and a 5 (q 1 ,q 2 )-f a^fa^qg) for all q ]L ,0 £ q 1 6 5, and all q^O ^ q^ 9. 
A formulation, neglecting input-output symbols, is as follows! 



V0q L 5 



b = f b,b„, 
o ol2' 



b l= + a i*i a 2V 

▼ <V b 2 = + a 3W*«iV * 
The quantification of q, in the master phase equation will be set up first. 
Then h, will be programmed, since it precedes b„ in the phase equation. 
Before the program for b g is composed, the quantification of q» is set up. 

Remark ! Each quantification corresponds to a 'loop' in the flow chart. 
In example 8.3, the loop for q^ is traversed six times. Within the q, loop 
is the loop for q^, which is traversed ten times for each traverse or the 
q, loop. Six values of b. are computed in the q, loop, whereas sixty values 
of b 2 are computed in the q^ loop. 

e.g. 8.4 . It is required to compute a-(q- )+» apCcu) and sinfa^q )^~ a^q,)) 

for all q-,0 4 q, ^ 5. A formulation, neglecting input-output symbols, is 
as follows! 

V'.V 

VOq.^ as Sin b 2 , 

b 2 X 4 * a l q -l a 2%' ' 

Here, the master phase equation serves only as a starting point. Since 
the b, equation is quantified, all dependent variables in that equation 
are also quantified by that same quantifier. Thus, the bp-equation is 
implicitly quantified. 

As a rule of thumb for specifying the order of computation, the 
formulator can use the following brief outline of the operation of the 
Encoder. The Encoder will start at the master phase equation, scan right 
for the first dependent variable, b- in example 8.%. It finds the b, -equation, 
sets up the quantifier for q,, then scans the right side for b*s. If one 
is present, the Encoder finds its defining equation. It repeats this 
process until an equation is found which contains no b's on the right side. 
This equation is programmed first. Then the path of search is retraced, 
ignoring b*s which have already been programmed. In this example, after 
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the quantification of q, is set up, b will be programmed and then K . 
(Note that the quantification could have been written in the master j5ha.se 
equation.) 

Rule 8»^ » No quantifier can be written more than once in a formulation. 

Thus, not only is it unnecessary to quantify b in the preceding 
example, but it is an error, since this would require writing , VOq,5* twice. 

9. Branching 

We now introduce a logical operation into the language. This will 
permit the definition of a quantity that is to be computed by one of two 
or more alternative formulas, depending on one or more conditions. This 
situation corresponds to a branch point in a flow chart of a calculation. 

Def. 9.1. A branch equation is an equation of the form, 
ys P^Xg, *, O, , 

where y denotes the dependent variable being defined, p. is a special 
punctuation symbol called a •branch symbol * \ <| denotes one of three 
library functions for the conditions £ , < , and sr ", x, and x are variables; 

ij denotes the b-term which defines y if the condition <Jx.x 2 holds, and £ 

denotes the b-term which defines y if the condition does not hold. 

Def, 9.2. The definition 6.k of b-equation is hereby extended to 
include branch equations. 

Def. 9.3 . A branch r-equation is an equation of the form, 

where i denotes the dependent index being defined^ P^ is the branch symbolj 
(J» is as in def. 9.IJ x. and x are indices or independent variables I 
\|r denotes the index term which defines i when (Jx,^ holds, and £■ denotes 
the index term which defines i otherwise. 

Def. 9.4 . The definition 6.2 of r-equation is hereby extended to 
include branch r-equations. 

e.g. 9.1 . The branch equation, 

\ = p 4 f 10 a 5 7 *' + a l a 2' "*l a 2' 

15 



NAVOHD Report 4209 

defines b, as equal to (a--|- O if a~< 1, and (a, - O if a~^ 1, 
i.e. f. denotes * < ' and is a library function. 

e.g. 9.2 . r x s fy'lAV + q 2 1 ' "^J 
This branch equation defines r. to be equal to (qg-4- 1) if q-, ^ a^, and 
equal to (qg - 5) if q x > « 2 > I.e. f^ is the library function for * j£ ». 

Observe that multiple branches can be set up by using a succession of 
branch equations in which one of the two alternative terms is f.,x, where x 
is a dependent variable (or index, as the case may be) which is defined by 
another branch equation. 

e.g.9.3. bjS p 4 f io a 5 7,+ "lV f lV 
^2-- P l(.^lo a 5^' "*i a 2' * a i a 2* 
These equations define b^*. (a^-f ag) for a_< 7, b^s (a^ - a^) for 7 ^ a~< 9> 
and b^s a 1 «a 2 for a ^ 9. 

Remark . We could dispense with the branch equation as such. Suppose 
<}(x) is a function which is equal to 1 if x & a and equal to if x > a. 
The equation, y = (}f ■+• (Jg, is equivalent to a branch equation for y, with 
f and g as the two alternative terms. The function (J is a primitive 
recursive function, and could be incorporated into the library as one of 
the given functions of the system. Nevertheless, the branch equation is 
included in the language for practical reasons. Many mathematicians are 
accustomed to that terminology, and it leads to more efficient programs. 

10. Vector Equations. 

A vector quantity is an n-tuple of scalar quantities J i.e. the 
components of the vector. In ABES II, we shall not provide any special 
letter to designate a vector, since all arithmetic will be carried out with 
the components of vectors. In fact, if the components are denoted by 
different b*s, no new syntax ,is required at all. However, it may be 
convenient to represent the components of a vector by an indexed b, such 
as b^q., where q. xz 0,1,..., n. In that case, several situations can arise, 
each requiring different phasing, quantification and output specifications. 

If the components, b-q,, are defined by a single formula, then an 
output specification may be required in the equation for b^ This is 
determined by rules given in section 13. 

If the components, b cu, are defined by different formulas involving 
different arguments, and Ir it is necessary to compute all components 
before the computation can proceed, then a »vector equation* should be used. 
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Def. 10.1 . A vector equation is a set of consecutive b-equations 
^fHtfje form, 

^ s 09 4)' 
y -09 *1» 



09 



4 n > > 



where y denotes the dependent variable which represents the vector and 6 
denotes the b-term which defines the componentr y(i), i s 0,l,...,n. 
Each component must then be called for by writing y with a suitable index. 

e.g. 10.1. Given the acceleration vector (cos t, sin t, t+1), 
compute the force vector, b., for a mass a. at time t s a . An ABES 
formulation follows. 1 ° 

b o sf oW 

Vo *i 2b i = ,a iVi> 

(\ "09 cos V 
b 2 - 0Q sin a Q , 

b 2 =1 09 "*" 1 a o' ) • 

How, it frequently happens that (k + 1) scalar quantities, y(i),0$ i^ k, 
are defined by k-f 1 different functions, g., of the same arguments, x.,...,x . 

In such cases, it is customary for the mathematician to write, 

y(i)-=-S ± (\fo,\), 0$ i $k, 

and then define the functions g. in separate equations. This is especially 
convenient when another vector, z(i), is defined by the same functions 
with different arguments J i.e. 

z(i) -%(«!*. -.»%), O^i^k. 

This type of formulation requires the introduction of an indexed function 
symbol, and a means of defining it. 

Def. 10.2 . An indexed function is an expression of the form, gi, 
where g denotes an auxiliary function and i denotes an independent index. 

e.g. f 51 q x ," f^qg; 

Def. IQ.g . A vector auxiliary equation is a set of consecutive 
auxiliary equations of the form, 
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s - 09 <t> 

• • 

« 

o • 

8 =09 ^ } » 

where g denotes an auxiliary function, and the <L denotec free terms which 
define the functions denoted hy the indexed function, gi, &i$k. 



Rule 10 .X . An indexed function, gi, can be used only in a b-equation, 
g must be the only function in that equation. A vector auxiliary equation 
must be written to define gi. The index, i, must be quantified somewhere 
in the formulation. If y denotes the dependent variable being defined by 
the indexed function, we write the b-equation, 

y« gi ^ ... x n , , 

where the x«s denote the arguments of gi. If it is necessary to refer to 
several values of y in the formulation, one writes y with a suitable index. 
In that case, the equation for y must contain a storage symbol as expla ine d 
in section 13. 

As a simple example, consider the problem, 

y i~ h^ ±h ^ ±)h °^ i - 1 
ZjS efoii),*^*)), 0< j 4l, 

s o~ sin ( c i"*" c 2^ 
g;L = cosfc^ c 2 ), . 

Letting ^at y, bgS z, q x = i, g^st i, f^S g, W& Xj^ 1 )* 
"2*1- at 2^ ± )» 8 5%*' X 5(J)» a Vfe= %(3)> ™ obtain the ABES formulation, 

< f 6o - 09 Sin + C l°2' 

f 60=09 COS + C l C 2'> 
b o s f o b l b 2> 

Voq i lb i- f 60Wl¥l' 

Vpqglbg s f^a^a^, 



11. Recursion 

If the ADES language is to allow for the formulation of primitive 
recursive functions, the syntax must include rules for definition by 
recursion. This is probably the most difficult part of the syntax. 
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In ADES II, we provide for the definition of dependent variables 
by several different kinds of recursion equations. The various types 
of recursion equations are distinguished by means of identifying sub- 
scripts attached to the equal signs. With the exception of a vector 
equation, the equal sign in any non-recursive equation has the subscript 
00, which we have agreed to omit throughout this report. In recursion 
equations, a two-digit subscript, n_n , is written explicitly with each 
equal sign. The following table lists the value of il and n , and the 
names of the corresponding recursions. 2 

Subscript Type of Recursion 

1 1 Simple scalar preceding-values 

2 1 Simple scalar course-of-values 
5 1 Simple vector preceding-values 
k 1 Simple vector course-of-values 
5 1 Simple vector one-row 

1 2 Double preceding-i row 

2 2 Double two-preceding-rows 

3 2 Double one-row 

k 2 Double course-of-values 

3 Triple special 

The names are intended to describe both the structure of the mathe- 
matical formulation and the kind of store instructions in the program 
produced by the Encoder. 

The second digit, n , indicates the number of independent indices 
Involved in the recursion. 



e.g. 11.1. For example, the Newton iteration algorithm for the 
square root of a. is a recursion on one index. It defines a scalar function, 
and the new value in an iteration depends only on the preceding value. 
Hence, it is a simple scalar preceding-values recursion ] its subscript is 
U. In ordinary mathematical notation, it might be written as, 



x(i + l)i(l/2)[x(i)-|-a ]L /x(i)] , 
x(0) S 3, 



where the recursion index, i, goes from zero to an upper bound, u, which 
depends on the accuracy desired. We shall rewrite this later as an ADES 
formulation. 
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e.g. 11.2. Che summation, b, sa £""«„ ^(q,,), also can be written 

as a recursion with subscript 11, namely, 

V*L + X) * a l (q l )+ V*L>' O^in-1, 
b x (0) S 0, 

e.g. 11 .3 « As an example of a double recursion with subscript 12, 
we take the recursion formula for the binomial coefficients. Denoting the 
coefficients by b (i,j), we write in conventional notation, 

^(i + 1, a + 1) = b x (i, 3 + 1) + bjU, j ), 

b^H-1,0) *1, for all 1< -1, 
b^O, j + 1)850, for all $ Z 0, . 

Def. 11.1. A simple scalar recursion equation is an equation of the 
form, 

VLiux s^, * , ..., + ffl , , 

where x denotes the dependent variable being defined, i denotes the recursion 
(independent) index, <} denotes the b-term which defines x(i+- 1) as a 
function of x(i), x(i - 1), ..., x(L), and other variables, and f Q , ..., t m , 

are the b-terms which define the initial values x(L), x(L+ 1 ),..., x(m), 
respectively. VLiu is the quantifier for i. It is the only quantifier 
permitted in the equation. 

It is understood that the x on the left side of the equal sign denotes 
x(i+ 1). If for each value of i, the term, d, involves only the m-V 1 
values of x preceding x(i+ 1), i.e. only x(i), x(i - l),...,x(i - m), then 
n, — 1, and we call this a preceding-values recursion. If <} involves any 

other set of values of x(e.g. all the values prior to x(it 1)), we call 
this a course— of-values recursion and write a. = 2, 

Remark I In a course-of-values recursion, the Encoder will reserve a 
total of N(l) computer storages for the values of x, where N(i) is the 
M vrtjmm number of values which the index i assumes, whereas in a preceding- 
values recursion, it allots only m-f2 storages. Thus, although 1^= 2 
will always yield a correct program, n,= 1 should be used whenever possible 
since it conserves computer storage. 
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e.g. 11.4. The Newton iteration algorithm of e.g. 11 .1 for Ta. sz *f bl. 
can be formulated in ADES as follows (neglecting Input-Output). 

o o 1' 
Vo qi 6 b x m n / +b l4l / ^W" f l 5 > ' 

Here, six iterations are specified. Later, we shall introduce a device 
to permit the number of iterations to depend on the accuracy desired. 

Def. 11.1a. The class of b-equations is hereby extended to include 
simple scalar recursion equations. 

A system of simultaneous recursion equations on one index will be 
called a * simple vector recursion 1 . 

e.g. 11.5. Suppose the quantities x and y are defined by the system 
of equations, 

x(0)35 a^ 
y(0) s a 2 , . 

The vector (x,y) is thereby defined by a 'simple vector recursion*. The 
formal definition is as follows. 



Def. 11.2. A simple vector recursion is a system of equations of the 
form, 

VLiu(y S„ 1 4 Q , ir^, V 01 , ...... V^, 



n,l 

yi-^1 <*!> . v * tt * ••» *W 



^=^1 A' *n0' *m' ••" V* > » 

where y,., ,.y_, ..., y denote dependent variables representing the components 
ox n 

of the vector. There are two possible cases* Either each y denotes a 
distinct dependent variable, or all the y's denote the same dependent 
variable. In the former case, the ith value of the jth component is 
denoted by y,i. In the latter case, it is denoted by yi^, where I j^ n, 

j 
and is defined by the jth equation. 

It is understood that y_, ..., y on the left side of the equations is 

u n 

an abbreviation for the values at the new point, i + 1. Thus, <L defines 
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y.(i + 1) in terms of the preceding values y Q (i)> ♦••> y^ 1 )* yo^ ~ ^ •••' 
y n (i - 1), ..., y Q (L), ..., y n (L), and the new values y k (i + 1), where k<. j. 

*. » .... ty* denote b—terms which define the initial values y.(L), .... 

jo' * jm 

y.(L+ m), respectively. As in a simple scalar recursion, «V Liu* is the 

J 
quantifier for the recursion index, i. 

If each (f . involves only the m-f- 1 preceding values of y Q , ..., y n , 
and the new values of y fc , k< j, the subscript i^ss 3 in all equations. As 
a special case, if each (J . involves only y fc (i + 1) for k < J, and y k (i) for 
k > d, the subscript iuse 5. This is the vector 'one-row* recursion . In 
all other cases, il s ij i.e. they are treated as course-of-values vector 
recursions . 

No quantifiers are permitted, within the parentheses. 

In calling for a vector recursion, one refers to the desired component. 
Reference to any component causes the entire recursion to be computed. The 
components are computed in the order in which they are written. 

e.g. 11.6. We rewrite example 11.5 as an AKES formulation. 

b «■ f b^ , 
Voq^ = ?1 + VAV f lV 

b 2=31 /b l 4 l b 2V f lV > • 

In the above formulation, the vector is denoted by (b x , bg). 

The problem can also be formulated with the notation (b^O), b^O-)) 
for the vector. Thus, 

r l = ' f l " r 2 sSf l 1 -' 

b o ~ f oV 
YOq^O^ « 31 + b 1 4 L P 1 W 2 , t A , 

b i= 5 i / ¥i r i¥iV f iV > • 

In this second case, the components involve two indices in the recursion 
fo rmulas . The final results, however, when referred to in another part of 
the problem, are denoted by hji^ and b^. 

If the components of the vector in a recursion are denoted by hj(j), 
£ i & n, and if the same formula with different operands defines atL 
components, then the simple vector recursion can be regarded as a double 
recursion. Thus, using i as the recursion index, b L (i+ 1, J+ 1) would be 
defined in terms of the preceding values of b^ The final components of 
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the vector would then appear as the final line of the double recursion, 
and would have to be referred to as such. This will be further explained 
in the ensuing paragraphs on double recursion. 

As shown in example 11 .3, a double recursion involves two indices, 
i and j. We shall speak of the (i,j) , grid* and' refer to •points* in this 
grid in an obvious way. A double recursion defines a quantity, b. say, at 
the grid point (i "1-1, ,3 4-1) in terms of its values at grid poinds which 
•precede* (i + 1, j -hi). Following Peter £33, we say that the grid point 
(a,, P x ) precedes (a 2 , P 2 ) if either a 1 < a 2 , or 0^= a 2 and p x < p 2 « 

Geometrically speaking, grid point P, precedes P g if either P-j^ is in a 

lower row than P or if P.. is in the same row and in a column to the left 

of P p . i and j will be called the row and column index, respectively. 

The order of computation in a double recursion will be along a row from 
column to column to the last point in the row, then up to the first point 
in the next row and so on. 



Def . 11.3. A double recursion is a system of b-equations of the form, 



V Li W y l = 1^2 <*1 

xl.: 



y 2 ~ n, 2 %' 



m 9 

y k-l "^ 1^2 1k-l' 
^ L 2 i 2 U 2[ y k r k 1 = : n 1 2 ( S ! k> 

■'..•■• 
•''■■• 

Vk t «^2 <*s'] 

y s + 1 ==^2 4 S + 1' 

y m m r^2. ^m* ) • ■. 

i_ and i p denote the two recursion (independent) indices. i 1 is the 
row index. i p is the column index. (The column index, i 2 , will run from 
its lower bound L p to its upper bound u for each value of the row index, 1^) 
y*y a.1* •••» y denote dependent variables which are to be computed for 
all i, and for all ± p . The equations for these variables are enclosed in 
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the square brackets, y x , ..., y^ and y fl ^ ..., y m denote b«s which 

are to be computed only for all i,. Tbe equations for these b's are 
within the parentheses, but outside the brackets. 

For k £ j £ s, in the equation, y. s 2 <L, <L denotes the complete 

recursion formula for computing y. at the grid point (i,+* 1, i 2 +" *)• 
Unlike simple recursion, the equation for y . will usually be a branch equation 
in which the alternative terms define the initial values y,(0, !„■+• 1) and 
y-(l 1 + ijO)* «« d *ne general value yA.\+ 1> * 2 + 1), . Because of the 
choice of the indices i^ 1 and !„+ 1 for the •new* grid point, the 
lower and upper bounds in the quantifiers of i. and i must be one less 
than the actual bounds of i. and i J e.g. if ^ i. £ a. + 1, we must 
write • V-li^, » to insure that y,(jL + 1,0) starts with y 1 (0,0). 

Of course, for j < k or j > s, there is only one pertinent index, 
namely i^. Hence, (J defines the initial value y.(o) and the value y (i.-^. 1), 

again by using a branch equation. 

Note^ Wo quantifiers other than those for i_ and i are permitted 
within the parentheses. e ~ 

Def. 11.5 is continued in the following rules, 11.1 to 11.5. 

Rule 11.1. The equations in a double recursion will be computed in 

the order in which they are written J i.e. y. is computed after y. _ and 

j «}""•>- 

and before y. . n . This means that 6 . can contain operands of the following 
types: J+1 J 

(1) y m (a>P)> where a < J, and k £ m £ sj a £ 1^ + l,p £ i 2 + 1 . 

(2) y (a,p), where m,2£ j and k £ m £ sj and (a,p) precedes 
(i x + 1, i 2 + 1), ^fe 0, i 2 2.0. 

(5) y (a), where m< j, and m<korm>sJa£i,«* , l 

(k) y (a), where m £ j and m < k or m > s, oc£ i~, (i, > 0). 

Rule 11.2. In a double recursion, the digit, a,, in the subscript of 
the equal sign of the equation for y . is determinedly the row indices 
that go with the occurrences of y. in the (J*s. Each equation may have a 
different n subscript. J 
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(1) If for all m < j, ^ Involves only y.Ci^p)/ where 3 is 
subject to rule 11.1, and if for all m> J, <L Involves only y A (U 4" 1,P), 
where is subject to rule 11.1 and if q involves only y (1 4- 1,0), 
where p < i & 4* 1, or ^(i^P)* vhere p** ig -V 1, then x^s 3 in the 
equation for y.. 

(2) If only 1. and i.«f* 1 appear as row indices of y,, then we 
write n.s 1 if we cannot write n, s 3. 

(3) If only 1,-1, i^ and 1^4- 1 appear as row indices of y>, 
then we write n. ss 2. 

(It) In all the remaining eases, we write n.8 4, 

The above rule for determining n. will assure an economy of storage 
in the computer. However, if economy is not an Important factor, the 
formulator can simply write a, S k in all equations. Each equation will 
then be treated as a course-oxUralues recursion. 

Role 11.3. If n.ss 4 in an equation within the brackets, the formulator 
must write a dependent index immediately "to the right of the dependent 
variable on the left side of the equation. The precise nature of this 
index will be explained in Rule 13.6 in section 13, on Input-Output. Its 
purpose is to provide store instructions. 

Note I Role 11.1 regarding the order of the b-equatlcns within the 
parentheses of a double recursion must be strictly observed. This rule 
prescribes necessary conditions on the order, but it does not determine 
the order uniquely in all eases. Hence, some of the ordering is at the 
formulator *s option. 

Rule 11.4. If one of the equations within the parentheses involves 
as an operand a dependent variable, b., which is not part of the simulta- 
neous double recursion, then the equation for b - should be written outside 
the parentheses. 

As remarked earlier, each quantification corresponds to a loop in 
the flow chart. Thus, in a double recursion, the quantification of i, 
sets up a recursive loop which is repeated for -Lsi-S u,. Within 
the 1, -loop is the recursive loop for iL. Recursive*!©©!* differ from 
ordinSry loops in that the quantities computed are stored in a recursive 
manner. Thus, recursive loops allow several interdependent variables to 
be computed simultaneously; i.e. if b is a function of b and b_ is a 
function of b_, this Implies that a recursion is taking place, since 
ADES does not permit implicit functions. 
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To provide greater flexibility in the ordering pf recursive loops, 
the double recursion definition is extended as follows. 

Def. 11.4 . A general double recursion is a double recursion in 
which there is more than one column index. Each column index is quantified 
separately, and the b-equations to which it applies are enclosed in square 
brackets . 

For example, a general double recursion in which there are two column 
indices would have the following format I 

* o 

* • 

y k-l s il2 Tfe-J.' 

• » 

a . ■ • 

y s ♦ l m ^2 ^s +■ 1* 

y m *n-2 m' 
V -u 5 u 5 [y m + ^a ^2^+1* 

y h+ 1* 1^2 ^h+ 1' 



y. 



t ~n,2 



tv ). 



Rules 11.1, 11.2, and 11.3 en order and indices apply, with obvious 

modifications to account for 1,. 

3 

e.g. 11.7 . we reformulate example 11.3 In ABES (omitting input-output 
specifications ) . The binomial coefficients are denoted by b_, the row 
index by q_, and the column index by q^. 
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r l~ "*" < 32 1 * 

V-lq. 1 a 5 (V-lq 2 a 5 C.\ ss 12 P^ 1 o a 2 ' -'l 1 " *lVl * 

b 2~ ^10*1°' f l°' + ¥l r lVlV • 

b o* f o b l' 

Remarks J There is only one recursion equation, namely, that which 
defines b, . Since there are two initial conditions, a triple branch is 
required. The first branch condition, f_ a 2 0,is *&„< 1 , where &^zz q~. 

If a < 0, i.e. if q^S — 1, the equation defines b., to be equal to 1. If 

a o £, 0> n, ss b_, where b_ is defined by a branch equation outside the 

parentheses. The branch condition for b 2 is , a.< 0*, where a.s q,. If 

a, < 0, i.e. q 1 s — 1, then b-S b_— 0. If a,2 0* we have the main recursion 

formula. Part 2 of Rule 11.2 applies, and we write 12 as the subscript of 
the recursion equal sign. Therefore, this will be treated by the Encoder as 
a preceding-row recursion. (a + 2) storages will be reserved for the 
•preceding* row corresponding to the index q_, and (a + 2) storages for 

the •current 1 row corresponding to q.-^ 1» Note that in applying Rule 11.2, 
we must check through all b-equations referred to by the recursion equation 
for b 1 . 

e.g. 11.8 . In the Choleski method for solving systems of linear 
algebraic equations, a square symmetric matrix a,(i,j), of order (n+1) 
is given, and a triangular matrix b_(i,j) is computed by the formulas, 

i > J, b^i,,})* £a 5 (i,j) - 4* 1 O b l< i ' k > b l ( ' 3 ' k i] Aam>> 

is,} b 1 (i,i) = ^"[a 5 (i,i) - 4"io b l (i ' k > b l (i ' k '* 

Denoting i by q x , 3 by qg, (n-1) by a^, k by q^, the expression in 

brackets by b , and the summation by b , we can formulate this in ABES as 
follows (again without input-output)! 

'3.C+4L 1 * r 2- + ^2 1 ' - 

V-lq^C V-lq 2 q 1 [b^ ^VloVV/VlVV W J ) 

b 2 S " a 5 r i r 2 b 3' 



V^V^ii+^iV^ViVW f i°" 



b as f b. , 
o o !'■ 



Remarks! m the input specifications, a^ and a 2 are identified as 
q and q_ respectively, i.e. they are integer-valued variables (Def.6.5). 
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i 
This artifice must be used since a b-equation cannot contain indices 
as operands. Since all values of b. are involved in the recursion, it 
is considered to be a course-of-values recursion. Hence, the subscript 
is te. The symbol r, is a storage (i.e. output) specification. It will 
be discussed in section 13, Rule 15.6. 

In the preceding section on vectors, it was said that a simple vector 
recursion in which the components of the vector are defined by the same 
function can be formulated as a double recursion. The next example 
illustrates this. It is the well-known Gauss-Seidel method for solving a 
system of linear equations. 

e.g. 11.9. The system of equations, 

s'L^Ci,^) b^j) = a 2 (i), is 0,1,. ..,n, is to be solved for the 

vector, b-Cj), given the matrix a^ijj) and the vector a 2 (i). Denoting 

the iteration index by k, the Gauss-Seidel formula, in conventional 
notation, can be written as, 

- E q = 0*1^ + ^V* 4 * ^aj 

This is a vector recursion on the index It. It is of the type in which all 
components of the vector are defined by a single formula. As mentioned 
in section 10, such a vector recursion can be treated as a double recursion. 
The Gauss-Seidel formula illustrates this. In ADES language, if we write 
q for k, q^ for j, q for q in the first sum, which we denote by b g , 

write q. for q in the second sum, which we denote by b,, and a, for n, we 

get the formulation (using a^4* 2 iterations)! 

o o 1* 
r l s * 1,q 2' r 2 = + lq l' r k S + 2 'V r 5 = "V"' 
V -Iq^ (V-l V 3 L b l *12 / "Vl+Wl'l ']') 
V r ^% a 5 b 2 S 11+ b 2% B VlWlV f l° ' 
Voq^qg ^ = 11 +^ # YM r 2V f l ' « 

Next, we consider the following problem. Suppose that the elements 
of a two-dimensional matrix, ^(0^,0^), are to be computed by a simple 

recursion on q. . In some cases, one could set up the quantification to 
take place in the order q , q^q,, and a simple scalar recursion formula- 
tion would suffice. Thus7 for each element, b^g^o^), a simple recursion 
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would be performed. Even if this procedure were possible mathematically, 
it might not yield an efficient program since same quantities might be 
recomputed many times. In cases where the recursion is simultaneous, 
that is, b 1 (q 1 + 'i*9gMj) depends on some b^q^ocp) where (oc,0) # (q^qj, 

the simple recursion procedure is no longer valid. The quantification 
must now take place in the order qgjq^q, or q-jqu^q, > that is, for each 

q 1 , we compute b^ for all q* and q^. This is, in effect, a triple recursion. 

In ADES II, a special type of triple recursion is provided. It will 
cope with the above matrix recursion, and with any triple recursion which 
satisfies the special requirements explained below. 

Def. ll.g. Let the symbol "GBR* denote a set of b-equations written 
in the format of a general double recursion, except that all the equal signs 
are written with subscript 03 • A special triple recursion is a set of 
b-equations of the form, 



v^iM/i-c^ 



y k S 03^k' 

( GBR )^ 

y k4 1 *03 ^k-*- 1' 



ys-ojV 

(GDR ) 2 

y m S 05W J ' 

Thus, a triple recursion consists Of b-equations interspersed with GDR*s, 
all quantified by the i_ quantifier. In each GDR there is one row 
quantification and any number of column quantifications (possibly none). 

We shall not attempt to give formal rules for the writing of indices 
in a triple recursion. Instead, we shall explain briefly the order of 
computation and the manner of storage in the computer. This will indicate 
how to index the dependent variables. 

Let i denote a row index which precedes the parenthesis in a GDR. 
Let i be a column index which precedes a bracket in that GDR. If y, is 
defined by an equation within the i, brackets, the recursion for y. proceeds 
along the grid points of the (*-,+ *■) -plane, running along the (ip+ 1) -row, 
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as the column index, i,, runs from its lower bound to its upper "bound. 
The storage of y. values is limited to at most a two-dimensional array . 
All storage specification is accomplished by writing a suitable 
dependent index to the right of y . in the y . -equation. This is 
explained in section 13. If y . is stored as a two-dimensional array, 
this implies that when the index i. increases and a new value 
y (i +• 1, a,0) is computed, it Immediately replaces y 1 (i l ,o,3) In 

storage. With this convention, only the row and column indices need be 
specified when y . is used as an operand, y . can also be stored as a 

one-dimensional array if the problem permits. In that case, only one 
row of y. is stored in the computer at any juncture. When i g increases, 

the new value y.(ip+" 1,0) replaces y.(ig,p) immediately. Hence, if y, 

is used as an operand, only one index need be written in this instance. 

If the y, equation occurs within the ip parentheses of a GBR but not 

within any of the brackets, the recursion for y will run from row to row 
in the (i,+ 1) -plane as i 2 varies. Here, we Have the possibility of 
storing y* as two-dimensional (i.,i 2 ) array, or as a one-dimensional i g 

array. In the latter case, when 7*1*^+ !*«) is computed, it immediately 

replaces y.Ci^a) in storage. Thus, only the array y*(ig) is in storage 

at any time. Part of this array may belong to the (ij+ 1 ) -plane and the 
rest to the previous plane. When y, is used as an operand, only one 
index is required in this case. 

Finally, if the y, -equation occurs outside all GDR*s, it is computed 
only for all i, . Here^ the new value of y. does not immediately replace 
the previous value. This replacement is done after all y*s have been 
computed in the i. loop. Thus, after y. is evaluated in the (1^+ l) -plane, 

there are two values of y. in storage, namely, Y*(i L +' *-) and y^ 1 *)- *° 

use either as an operand in a formula one simply writes the proper index, 
keeping in mind that the y«s are computed in the order in which their 
equations are written. If all values, y (j^), are to be stored, an index 
must be written as explained in section 13 . 

No quantifications are written within the braces other than those in 
the GER's. However, if it necessary to introduce other quantifiers, one 
can write a phase equation, y s qx^^i' at the appropriate place within 

the braces, and then quantify the y. -equation, which is written outside 
the braces. Thus, for example, y. might itself be defined by a triple 
recursion. The storage for y ± is specified in the y ± -equation. 
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12. Minimization. 

In section 3, brief mention was made of the minimization operator 
which is introduced in the theory of recursive functions to eliminate 
multiple recursions. We incorporate a form of this operator into ABES II 
as a library function. However, we do not propose to eliminate multiple 
recursions. Rather, minimization will be used frequently within recursion 
formulations to define the upper bound of a recursion index. It will also 
be used in 'table look-up* operations. 

Def . 12.1. A minimization equation is a special equation of the type 
i x f^ixy, , 

where j is a dependent index, i is an independent index, x is an indexed 
independent variable or numerical constant, y is a dependent variable, and 
f u denotes a special library function known as the 'minimization operator* . 
Tne minimization equation can be read as « j is the minimum value of i,i £ x, 
such that y £ 0' • 

Def. 12.2. The class of r-equations is hereby extended to include 
Trrtriimlzfl-tion equations . 

e.g. 12.1. In example 11.4, the Newton iteration algorithm for the 
square root was formulated in ADES language for six iterations. We now _^ 
impose the requirement that the error in the square root be less than 10 , 
if possible, but that at most six iterations be used. The formulation, with 
a minimization equation, becomes 

r i=Vi 6 *V 



r 



2 



= 4- qjl, 



b s fh^ 

b A S — f . -b_r o b n q.l0 , , 
2 abs. 1 2 r*! ' ' 



where we have written f . to denote the library function for the absolute 

_ aDS. 

value. 

Remark: The above use of the minimization operator in scalar recursions 

is so common that some special abbreviation would be feasible. We might 

simple write, . 

r, s f 6, x 10° l,x 10 , 
1 H 

it being understood that the first number, 6 x 10°, is the upper bound of 

the recursion index, while the second is the tolerance for the absolute 

difference between successive values of the quantity being computed. In 

vector recursions, the tolerance might be applied to the sum of squares 

of differences of successive components. 

31 



NAVDRD Report ^209 

The *table look-Hip 1 aspect of minimization is illustrated by the 
next example. 

e.g. 12.2. Let fifty values, a,(q-), of a function be tabulated 
against the corresponding values of its independent variable, a,(q,).. 
It is required to evaluate the function at some point, &^ t by 
linear interpolation. In conventional language, this value is given by, 

a ^ r i^~ x > " a ^( r i) r T 

\* a 3 (r l^ ^- 1) - a^(r') U " ^*1>J ' 

where r is the minimum value of qu such that ^.(q,,) -a p io. In ADES, 
the formulation (except for Input-Output specifications), is as follows! 

b Q s fb^ 
Other examples will be found in Appendix III. 



15. Input and Output Formulation . 

The mathematical and logical part of the formulation of a problem 
in ADES II, as described in the preceding sections, consists of r— equations, 
b-equations, and auxiliary equations. To complete the description of the 
language, syntactical rules will now be given for the specification of the 
input of data, and the output of results. 

Let a , a., a 9 , ..., a be the independent variables which occur in 

the r-equations and b-equations of a problem. The a^. represents; the 
different types of input data. Let H fc be the maximum number of numerical 
quantities to be supplied in the computer as data for a. . In problems 
which are toobe repeated for j different cases, the actual number, N, ., 

of data for a. may vary from case to case. Thus, N. ~ max (if. .\ . 

Def. l?.!. In what follows, we shall say that a variable, a^or b^), 

»is of degree two in a formulation* if it occurs at least once as part of 
an indexed variable of degree two; e.g. a^q^or \%i^U * f a variable 

a. (or b. ) is not of degree two, and if it occurs at least once as part of 

an indexed variable of degree one, then we say that a. (or b^) «is of 
degree one 1 . Otherwise, a. (or b. ) *is of degree zero*. 
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RemarkJ If an independent variable is of degree zero, one, or two, 
then it always occurs with zero, one, or two indices, respectively. The 
phrase 'at least once* really applies to dependent variables, since they 
can occur with a different number of indices at different places in the 
same problem. 

For an a. of degree zero, N k = lj i.e. one numerical datum will be 
loaded into the computer for a. . (Exception! If a. is an integer-valued 
variable (Def.6.5), no data is supplied.) If a. is of degree one, then 
N. > lj A. •column* of numerical data, to be substituted for a(j), Ofijiu, 
will be loaded into the computer. * 

If a^, is of degree two, a matrix of data, to be substituted for a v (i,j), 
will be loaded into the computer. * 

The loading of data into the computer* s internal storage in ADES II 
is under control of a 'Read* program composed by the Encoder. The Read 
program is executed by the computer just after the entire program has been 
loaded into internal storage. The Read program transfers data from external 
storage (e.g. tape, cards) to predetermined internal storages. In ADES II, 
the Encoder composes the Read program on the assumption that the data in 
external storage is arranged according to the following conventions. 

Computer Convention 1 . All the data for each a, is in a block of 
consecutive external storages. The beginning of theHjlock must be identified 
as belonging to the variable a, . If tape is used, the end of the block 
should also be marked. If punched cards are used, all data on a card should 
be for the same a. , and should be in the order dictated by conventions 2, 3, 
and h below. Furxher, all the data for an a. should be on consecutive cards, 
in order, and suitably identified. Just wha£ constitutes suitable identi- 
fication will depend on the particular computer. 

Remark! It is not difficult to see how these conventions for the 
computer can be rela x ed to allow for variations in the arrangement of input 
data. The seheme suggested here somewhat simplifies the Encoder and the 
Read program. For those problems in which the amount of data is somparatively 
small, the conventions impose little or no inconvenience. 

Computer Convention 2 . If a. is of degree zero, only one datum is 
supplied, i.e. this variable requires a block of external storage of length 1. 
The block should be suitably identified at beginning and end. The blocks for 
all a. of degree zero should themselves be consecutive. 

Computer Convention 3 . Following the data for all the variables of 
degree zero is the block of data for the first variable of degree one. 
All the data for a variable a. of degree one should be loaded in a block 
of consecutive storages in an order corresponding to a. (0), a. (l ),..., 
(ST — 1). Following the block of data for the first variable of degree 
e Is the block of data for the second variable of degree one, etc. 
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Computer Convention h Following the block of data for the last 
variable of degree one is the block of data for the first variable of 
degree two, a. . This is actually a matrix of data to be substituted for 
a.(i, j), < I < u , L ? ^ j i=- u . The linear ordering of the data in 

this matrix into a block of consecutive storages must correspond to 
\(0,L 2 ), aJOjLg-*- 1),..., £^(0,^), 8^(1,1^), a^l',1^-4- 1),..., a^l,^), 

..., & \S u x> L o^ ••"' \^ U 1 ,U 2^° In short > the raatrix should be loaded 
into external storage by consecutive rows. This ordering can be described 
as a mapping of the two-dimensional (i,4) grid onto a one-dimensional r-grid 
extending from some point a to the point a 4* N fe — 1. For example, if 
0^165 and < j < 5, the square 6x6 grid is mapped into the linear 
grid by the mapping, 

r = a +(j + 61), . 

Again, if ^ i i 5, i £ j 5 5, the triangular grid is mapped into the 
linear one by, 

If the matrix structure is more complicated, the mapping may be defined by 
using a tabular function. The data for the other variables of degree two 
is loaded similarly, and in consecutive blocks of storage, each suitably 
marked. 

The structure of the matrix of data for a variable, a. , of degree 
two determines a mapping according to Convention k. The formulator must 
write this mapping in the formulation as an index equation defining a 
dependent index. Such a dependent index will be called a *storage index* 
of type A, since it describes how the data for a is to be arranged in 
storage. However, in writing the , storage index equation*, the formulator 
need make no reference to any storage address J or equivalently, the initial 
storage address, a, is always taken to be zero. The storage index is 
defined as a purely mathematical function of two independent indices I 
q_o, representing the row index, and q^, the column index of a matrix. 
Tnus, for the square matrix with row index £ q^ ^. 5 and column index 
£ q^Q £ 5, the storage index equation is, ^^ 

r x = (qug + 6qo8>, (common notation) 

r x = 4. q« 6.9^, • (ABES notation) 

For the triangular matrix with 0£ q-g < 5, q^g < 3^ 1 5, the storage 
index is, 

»1 = % 9 + & - (6 = **ffi = W > , (common) 

r l - +q 99 " 21 " t* -6 * q 98 " 7 *%8 2 ' ' ( ATJES) 
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For the triangular mati^, O^q^^ the storage index 

equation is, 

r l = % + %8 (q 98 + X >/ 2 ' (common) 

r l r + %/ * ^8 + ^8 lo2 * ^^J * 

The formulator must also list the type A storage indices in a table, 
called the •Computer Table', which is part of the formulation. The 
Computer Table will be defined precisely in Def. 13. <f. 

Since the data for most of the independent variables of degree two 
in a given problem will usually have the same matrix structure, the same 
storage index will apply to most of the a*s. Rectangular and Triangular 
structures are the most common. Several rectangular arrays of different 
dimensions may occur in the same problem. In that case, it is more 
efficient to have a single function, f say, to denote the formula for 
the rectangular mapping. If f„. is a library function, the formulator need 
only write an equation of the type, 

r l s t r q ^ fi 99 a l* > 
where q-o and q^ are the row and column indices, respectively, and 
£ q-Q £ a, . If f is not a library function, the auxiliary equation, 

f r =+C 2' C lV > 

must also be written, but this auxiliary equation will serve for all other 
rectangular storage index equations. 

It should be evident that no particular row and column indices can be 
specified in a storage index equation for an independent variable of 
degree two. It is only necessary to distinguish the row index from the 
column index. With all the above considerations in mind, we adopt the 
following definition and rule for ADES II. 

Def. 1%2. A storage index equation of type A is an r-equation of the 
form, 

3- ( i%8 q 99 X 5 °" x n r » 

where J is a dependent index (called a type A storage index ), <\ is a 
function symbol, g-o always denotes the row index, q^ the column index, 

and the x»s are indexed independent variables. The function (J defines a 
mapping of the grid points (q^,^) onto a linear grid. If (f is not a 

library function it must be defined by an auxiliary equation in which the 
free variables \> c 2 f'> c n are *°' be Identified ^^ ( 3oB" 3 99 ,X 2 , *** ,X n , 
respectively, g^g and q™ are never quantified and can be used only in 
storage index equations Of type A. 
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Rule I?»l« To each independent variable, a. , of degree two In a 
formulation, there is assigned a storage index ox type A. A storage 
index equation of type A most be written to define this storage index. 
The function 6 in this equation depends on the structure of the matrix 
of data for a. . (J mus * ^P * he 'S^id points* of the matrix into a linear 
sequence according to computer convention k. The same storage index can 
be used for different variables if the matrix structures are identical. 
The type A storage indices must be listed in the 'Computer Table 1 ( see 
Def. 13 .4 below). 

Storage index equations are also written to specify the storage 
structure of computed results. The rules for these equations are somewhat 
lengthier. 

Def. 15.5 . A storage index equation of type B is an r-equation of 

the form 



jsf,i 



lV 



or of the form 



j — (fljigX^ • . .x n , , 

where j is a dependent index (a storage index of type B), A is a function, 
i and i are indices, and X~, ..., x n are indexed independent variables. 

If (J is an auxiliary function, the free variables c.^,^, ..., c q in the 

auxiliary equation are to be identified with i ,i , ...,x , respectively. 
A branch symbol is not permitted in a storage Index equation. 

After writing the b-equations, the formulator ascertains the degree 
of each b-symbol(Def . 13 .1) by inspecting the righthand sides of all the 
equations . 

Rule 13.2. If b is of degree zero, no storage index equation is 
required, b, will automatically be stored in the same storage each time 
it is computed. 

Rule 13.3. If b^ is of degree one or two, and if b^ is not defined 
by. a recursion or a vector equation , then a type B storage index, r * say, 
is assigned to t^. The symbol r . is written to the left of b^. on the 

left side of the equation for b. . A storage index equation of type B is 
added to the formulation to define r^. If \ is of degree one, the first 

form of the type B storage equation is used. The choice of the index i 
will cause h i , to be stored in separate storages for all i^ If b fe is 

of degree two, the second form of the type B storage index equation is 
used. The choice of the indices i x and ± 2 will cause b^ig to be stored 
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in separate storages for all i, and ig. The function (J is selected by 
the forraulator after he determines the structure of the (i ,i ) grid 
defined by the quantifiers of i_. and i , The discussion on tjifpe A 
storage equations is applicable if q^g and q__ are replaced by jL and 
i , respectively. "^y 

Note that the same storage index can be applied to different b»s 
if they have identical storage requirements . 

Careful consideration of Rules 13.2 and 13 .3 gives rise to several 
suggestions for writing b-equations. 

In the interests of economy in the use of computer storage, the 
formulator should avoid introducing unnecessary b-equations . m the 
computation of the wellr-formed formulas which occur as parts of a b— term, 
results are stored only when necessary . (This storing is automatically 
programmed by the Encoder. It does not concern the formulator.) However, 
the final result of a b— term is automatically stored whether or not this 
is necessary for the progress of the computation. For example, in the 
equation 

\zl • + a 1 a 2 + a 5 V 

the intermediate result, +■ a,a. , must and will be stored before the program 
proceeds. No other quantities will be stored, except the final result, b^. 
Now, if we reformulate this as follows, 

^ss+a^ag, 

b l s * b 2 b 3' > 
then h.,b p and b will be stored, b unnecessarily^ 

The indexing Of a dependent variable, h. should be avoided if possible, 
since this will require the Introduction of a siorage index. This will 
cause the Encoder to produce a program in which the different values of 
are stored in separate storages. A complete discussion of this point lee 
to rather involved and deep problems. We shall content ourselves, in the 
present report, with the formal rules and examples to illustrate what is 
at issue. _ 

5 L 

e.g. 13.1. Suppose it is required to compute -bj_« ^so '^(g^)* 

where b (q_) = ^i^) • %(%)« This should be formulated as follows. 

b = f o b x , 

^SS b l "tt+^ftV f l°' 
b 2 -• a l% a 2 1 1 , '• 
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Note that b„ is written without an index, since only one value of t> 2 at 
a time is needed, and the proper indexing is already contained in tne 
right side of the b_ -equation, 

a-1 
e.g. 13.2. Suppose b^s VsoV*lN^. + 1 ^ where a &^a. 

h ? (q.i) * shCq.!)* ^(li)* We shall formulate the problem in three ways, 
and mention a fourth. 
Formulation 1. 

r l " ~*J-> r 2 " **" lqL l' 

o o y 
Voo.^ b 3 SS U + b^ • a^a^a^a^, ^0, 

Formulation 2. 

b s f b-b,, 
w o o 2 2' 

^ 0q 2 a 5 r 3 b 2 s * *! V^' 
VOq 1 r^ 3 » 11 +b 3 q i .b 2 q 1 b 2 r 2 , fjO, 

Formulation 3. 

r l * "•j 1 * r 2 ~ + 1( ll' 
b Q a: fb , 

VOq^ ( b x S^h^, f x 0, 

*2 - 5 1*V 2 V2' • a i 0a 2°* 

*3 *51 +t 2 q l ,b l q iVl' f l ' } 

Another possible formulation would introduce a double recursion to 
compute b 2 . We shall not go into this. 

The preceding example is very informative and merits' study. In 
formulation 1, the symbol b 2 is dispensed with entirely. b ? is computed 
by a simple scalar recursion which Involves a rather long b&term. When 
this b-term is computed, no intermediate results will be stored. Thus, 
this formulation would be fairly economical storagewise. Timewise, however, 
it is obviously inefficient since each product of a^ by 8g is performed 
twice. 

In formulation 2, the duplicate computation of a»a 2 is avoided by 
computing all the products first and storing them for all qg. Then b, 
is computed by simply calling for the successive pairs of bg values. 
Since b_ is of degree one in this formulation, a B type storage index, r_, 
is written to the left of the bg -equation, and a B type equation, r^ «• Xj^t 
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indicates that b^q^ is to "be stored in separate storages for all q_. 
This formulation produces a faster program, but requires (a + 1) Inter- 
mediate storages. If a represents a large number, this formulation is 
uneconomical storagewise. 

In formulation 3, b is treated as the third component in a simple 
vector recursion. The second component is the product of a, by a . The 
first component is then defined merely as the previous value of the second 
component in the recursion, thus avoiding a duplicate calculation of the 
product a-»a_. The subscript 51 makes this a , one-row l recursion, which 
means that a minimum of storage is used^ i.e. the new value of a component 
replaces the previous value in storage. Therefore, this formulation is 
economical both timewise and storagewise. In fact, the Encoder will 
translate it into a program which is about as efficient as the program 
produced by a skilled human programmer, who will, in effect, program it 
as a vector recursion (although he may not think of it as such.) 

In formulation 2 above, the use of a B type storage index equation 
is illustrated. The equation is of the simple form, j » <fi , „ The next 
example illustrates the B-type storage equation with two indices. 

e.g. 13.3. Let it be required to compute three nxn matrices, 

b l^ q l ,q 2^ = V q l^ a 2^ f 

^(q^qg) s a 1 (q 1 )/a 2 (q 2 ), and their product 



b 5 (qyq^) = ^oV^CLj) W^' *** MES 



formulation is as follows (letting a, s» n-1). 

f 5 l-+V c l +lc 3' 



r l s V^V 



b o - W 5 ' 
V0q 1 a 5 V / 0q 2 a 5 b^ S f^, 

r 1 b 1 «+a 1 ^a 2 q 2 , 

Voq^Voq^a^ a f^, 

VOq 5 a 5 b 5 a il -hb 5 q 5 *b 1 q 5 q 5 b 2 q 5 Q 4 , ^0, 



Since b n and b« are of degree two in this formulation, the storage 
index, r , must be Introduced and defined as shown. It causes the Encoder 
to compose a program which will store each value of b^ and b„ in a separat 
storage according to Computer Convention 4. The program forb, will later 
call these values out as they are needed. Note the use of phase equations 
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and quantifiers to prescribe the order of computation. The phase equation 
for b must be introduced because the recursion equation for b 2 cannot be 
quantified more than once (Def . 11.1). 2 

Rules 12.2 and 13.3 apply to the storing of non-recursive dependent 
variables. We now state rules for recursively-defined variables. 

Rule 15.4. Let b fe be defined by a simple scalar recursion, b must 
be of degree one at least. If it is of degree one and occurs only with 
the recursion index, or with some index which depends on the recursion 
index, and with no other index, then no storage index need be assigned. 
' Sbe flnal v » lu e of b, in the recursion will automatically be stored like 
any dependent variable of degree zero. 

However, if h occurs with an index which is independent of the 
recursion index, then Rule 13.3 applies. The storage index pertains to 
the final value of b in the recursion. Since the recursion is repeated, 
there will be more tnan one final value. Each final value will be stored 
separately as prescribed by the storage index. The storage index is 
written immediately to the left of b fe , but to the right of the quantifier. 

If the formulator wishes to store one or more intermediate values of 
h. computed during the iterations, he must write the recursion as a course- 
or-values recursion. This automatically causes all values of b, to be 
stored. No storage index is required. * 

Rule 15.5 . If b, is one component in a simple vector recursion, and 
if the other components are not also denoted by b,, then Rule 13.4 applies 
to b k# If the components are denoted by b.(i), no storage index can be 
assigned. If it is necessary to store b,(l) for all j, this must be 
formulated by relabeling h^i), using the identity function. 

Rule 13.6. If b, is defined in a double recursion and its equation 
is outside the brackets, then Rule 13.4 applies. If the b, -equation 
occurs inside the brackets of a double recursion, then b. Is necessarily 
of degree two in the formulation. 

If the recursion is a course— of-values double recursion, then a 
B type storage index, r , is required as indicated in Rule 11.3. The 
index, r , is written immediately to the right of Tx on the left side 
of the equation, r is a function of the recursion indices ±.,± 9 » 
The formulator must define r by the mapping of the (i.,^) grid onto a 
linear grid, according to Convention 4. The structure or the (i.,,i 2 ) 
grid is determined by the quantifiers of i_ and i . Thus, r specifies 
the storing of all values of b (i_,i ). n 
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If b. occurs with one or two indices which are independent of the 
recursion indices, then Rule 13.4 applies for the storing of the final 
value of b • 

If b. is defined by an equation within the braces of a triple 
recursion, the degree of b. in the triple recursion should be determined. 
A suitable B-type storage index should then be written to the right of b. 
in the b, -equation. k 

Rule 13.7. If y denotes a variable of degree one which is defined by 
a vector equation (Def .10.1), then no storage index need be written. If 
y is a variable of degree one which is defined by an indexed function 
(Defs. 10.2, 10,3 and Rule 10.1), the y-equation should contain a B-type 
storage index. 

(Remark! The rules for storage indices, as given above, are to be 
regarded as temporary. It is planned to eliminate most or all storage 
indices in ABES IIIJ i.e. the Encoder of ABES III will take over this 
task.) 

Having written all necessary storage indices, the formulator can 
write the part of the formulation known as the Computer Table 1 . 

Def. 13 .4* Let x?, x°, ..., x? denote the independent variables of 
degree zero in a formulation, written in an order corresponding to the 
order of the data in external storage. (See computer conventions.) Let 

x,, x, ..., x , denote the independent variables of degree one. Let 
xi ', « v ..., x* ' be the independent variables of degree two, and 

(2) (o\ 

«J£ > • ••* 3 n *&e corresponding A— type storage indices. (The variables 

of degree one and two are in an order corresponding to the order of the 
blocks of data in external storage.) Let x. , x , ..., x denote the 

X c. S 

integers-valued variables, and i.,i , ..., i the respective independent 
indices with which they are to be identified (Def. 6.5.). Let j-,,.0.,3 
denote those independent indices which are used in recursions, and 
Mii )>••*> N (3 ) the respective numbers which specify the ma x i m um number 
of values each p index assumes. Finally, let K,,...,!^ denote the B-type 

storage indices and n(K-),...,N(K_) the respective numbers which give the 

maximum number of values of each index. The Computer Table is a table 
arranged as follows! 
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o 

*1 



o 



m 



,(2) 



1 
1 



4 



N_ 



m 

K (2) ,(2) 

* 1 ** 1 J 1 

c<!> h<S> *<!> 



Independent Variables 
of Degree Zero 



Independent Variables 
of Degree One 



Independent Variables 
of Degree Two 



(2) „(2) P) 

n n n 



i 
*1 

i 
*2 



n 



Serai— colon 



Integer-valued variables 



N(£) 



n(j p ) ; 



KJ. N^) 



Recursion Indices 



B— type Storage Indices 



*H N( V' 



where the various N*s denote actual numbers which specify the maximum 
number of pieces of data supplied for the adjacent variables. 



With the writing of the Computer Table, the formulation is virtually 
complete. It remains only to indicate the output^ i.e. which results are 
to be punched out on cards, transferred to tape, printed, etc. The details 
of output will vary with the computer. 
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In ADES, output is indicated by writing an output symbol, d, with 
suitable subscript, on the left side of the appropriate b-equation. The 
details of output are concealed in this subscript, and will not concern 
us here. It is sufficient to assume that a list of available modes of 
output and their corresponding subscripts is in the possession of the 
formulator^ e.g. d„ might mean •punch this result in the second field 
of a five— field caraj* d. might mean »print this result in the third 
field on the page and skip'' one line*, d^ p might mean *write and copy onto 
tape 2*, etc. 

Rule 13.8. If it is desired to have all the values of the dependent 

VvariaDle/Jb. J , as output, an appropriately subscripted output symbol, d , 

should be written to the left of lb J on the left side of the \b, 7 equation , 
If there is a B-type index to the left of b, , then d should be written 
immediately to the left of this index. 

In a recursion, this will cause only the final values of b to be 
processed for output. To effect output of all values in a recursion, d_ 
should be written to the right of b. and to the left of the equal sign. 

This completes the description of the input and output parts of a 
formulation* 

Ik . Summary. 

The Automatic Digital Encoding System is designed to enable the 
mathematician to present a mathematical problem to any modern electronic 
digital computer in a form closely resembling a Conventional mathematical 
statement of the problem. ADES eliminates the programmer, with his 
specialized knowledge of one computer. All the clerical tasks which the 
programmer would perform are taken over by a machine, the Encoder. The 
mathematician, or forjnulator, must first prepare the problem as a set of 
equations suitable for digital computation. He then formulates the problem 
in the ADES language. Tb±6 frequently requires little more than a trans- 
literation of the original equations. A complete formulation of the problem 
in ADES consists of the following; 

Computer Table, 
Auxiliary Equations, 
r — equations, 
b — equations . 

Itor the most part, the formulation of a problem can be planned and 
written with no specific computer in mind. It is true that the use of 
library functions does depend on the computer. However, some standard 
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assignment of subscripts for the elementary functions and other common 
library subroutines can be adopted \ e.g. f„ for + , f for — , etc. 
Failing that, the formulator can use the conventional symbols, and when 
a particular computer is selected, the conventional symbols can be 
replaced by the correct function symbols. The same procedure applies to 
output symbols, which are somewhat more difficult to standardize. 

When the formulation has been completed for a specific computer, it 
is punched in cards (or tape etc.) and loaded into the Encoder for that 
computer. The Encoder will translate the formulation into a program in 
the computer language, and will compile all necessary subroutines. The 
complete program with subroutines is then punched out on cards (or placed 
on tape) in a form ready for computation. This program can then be loaded 
into the computer. It will read the data into high-speed storage, perform 
the necessary computation, and yield the desired results as output. 
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APPENDIX I. Library Functions. 

Most of the common library subroutines can be referred to simply by 
writing the properly subscripted function symbols in the formulation. In 
system II, the computer library is assumed to comprise mainly these compara- 
tively elementary subroutines, and the preceding report is based on this 
assumption. Shis in no way implies that ADES is limited to such simple 
subroutines. As examples of more complex subroutines, we shall cite two 
situations which require somewhat special formulation. 

First, if the subroutine involves a recursion, and if there are 

operands containing the recursion index in the formulation, a separate 

recursion equation must be written. For example, suppose f denotes a 

subroutine for summation, Z. Then if we wish to formulate o?s Z a. fa ). 

1 q,ss Ol v i" 

it is clearly not sufficient to specify the quantity to be summed. One 
must also specify the index of summation (i.e. the recursion index) and its 
bounds. Further, if ZT _ Q^Cli ) occurs as part of another formula, it 

must be denoted by a b-symbol, which must then be defined by a b-equation. 
In this case, 



defines the summation. Note that this equation is not considered to be a 
recursion equation (i.e. subscript of equal sign is 00), since the subroutine 
will take care of all the recursive features except the quantification. 

The second situation which complicates a library subroutine can best be 
explained by an example. Suppose there is a subroutine for the numerical 
integration of any real elementary function, g(x), over a finite interval, 
£a,,a.]. The numerical formula would be of the type, 

where the a, (i) are constant coefficients independent of the function g, 

[aj-a,! 
and the limits a,,a^, while bg(i) depends on a, and a, . In fact,b_(i)=.a,W — sj-^-n 

Therefore, the integration subroutine depends on values of the function g 
-Men cannot be computed until b_(i) has been computed. The formula for g 
must be written with the argument, b„(i), and g itself must be identified 
by a dependent variable. A suggested procedure for referring to such 
library functions is the following. 

Let f,_ denote the subroutine which computes b.. sr J g(x)dx. To 
illustrate, let g(x) s sin x. The formulation for b. would be, 

Vc V \ * 'so VaVy 

b^ a sin b 2 , . 
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The above format can be standardized, that is, for each library function 
of the type described, the formulator writes a b-equation, quantifying it if 
an index is involved. On the right side, he writes the pertinent f-symbol, 
f , followed by indexed independent variables, which represent data to be 
supplied to the subroutine, and dependent variables . Of the dependent 
variables, the last one denotes the function to be operated on by f . The 
last b is defined by a b-equation. The other b-eymbols denote quan?Sties 
computed within the subroutine for use in the formula for the last b-symbol. 
They are not defined by b-equations . 

The Encoder must be slightly modified to cope with these special library 
functions. It oust recognize the f-eymbol before programming the last b. It 
must then operate so that the part of the subroutine which computes all b»s 
but the last is compiled into the program first. These b*s are replaced by 
addresses which are relative to the location of the subroutine in the program. 
Then the last b is programmed. The program for the library function is 
then completed. 
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APPENDIX II. 

The order in which the b-equations are programmed is determined by 
(1) the phase equations, (2) the parentheses in vector equations and 
recursions, and (5) by the order in which dependent variables are written 
on the right side of an equation. 

The Encoder determines the order of computation in a natural manner 
as follows. Starting with the master phase equation, it scans the right 
member from left to right for dependent variables. Let the first dependent 
variable which it finds be denoted by y_. Assume first that the y_ -equation 
is not a branch, vector, or recursion equation. The Encoder Scans the 
right side of the y -equation for dependent variables. If there are none, 
y. is programmed ana the Encoder returns to the phase equation to obtain 
tne next dependent variable. If the y, -equation contains a dependent 
variable, y p say, then the y„ -equation is inspected for dependent variables. 

If no dependent variables are found in the y 2 -equation, it is programmed 
and the Encoder returns to the y. -equation sad scans right again for 
dependent variables. If y 2 is the only dependent variable in the right 
member of the y_ -equationf then y. is programmed next, and the Encoder 
returns to the master phase equation. If the y -equation contains another 
dependent variable, y, say, the Encoder will proceed with y, as with the 
y -equation. Eventually, all dependent variables on the right side of the 
ylr -equation will have been programmed, and then y. itself is programmed. 

In this way, the Encoder traces its way through the formulation, 
programming the variables as they are encountered, whenever possible. This 
order is Interrupted, however, if a variable ,y. is one component in a 
vector equation. In that case, the Encoder will arrange to program the 
components in the order in which they are written, beginning at the left 
parenthesis and proceeding to the right. Similarly, in double and triple 
recursions, the equations are programmed in the order in which they are 
written. 

In a branch equation, the leftmost formula is programmed first. 
Normally, if a variable, y., is referred to in several formulas, it is 
programmed only once, at tie time of the first referral. An exception to 
this must be made in the case of a branch equation in which a dependent 
variable, y., is referred to In both branches. This is inefficient formu- 
lation, singe y. will have to be programmed twice. This situation can 
always be avoided by writing a vector equation in which y i is the first 
component and the branch equation defines the second component, for then 
y would be referred to prior to the branch and programmed once. Never- 
theless, the Encoder will accept the inefficient formulation, and program 
y. once in each branch. 
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APPENDIX 111 . Examples » . 

Each example will be stated first in a raore or less conventional 
mathematical way, and tlien a complete ADES formulation will be given. 

EXAMPLE 1 . 

Given the data points a, (q,), where o£ q-. 6: a 2 -4"l. 

It is required to compute 

+ 1 r , n% 



a i- 1 r -v 

b i=\ = .oti(%)J 



for all q_, where £ qu^ 2a . (The b_ are the elements of the normal 

matrix in the least square fitting of a polynomial of degree a over 
a_-V- 2 points,) Two ADES formulations are given. 

Formulation 1. 
Computer Table . %1 ^ ^ a^ r^gOO ^50 q^, 
Auxiliary Equations , f „ S. + ^'c^^ "+* c* 1 * 

r — Equations. r a "" r i 1 » 

r, = »2a , 
1 o* 

r 2 s f 50 r 3Vl' 

r^s -Vlqg, 
b '— Equations . b Q S f bib,, 

V I*. VfeW - i -V'oC b 2 r 2 **2 Vll a 5 - 1 ' ■ V> * ■ V 5 V3«2 1] > 
V°% r l.S*' f bV 

Remarks: ^(^i'Q^) '^{W^W' is defined by a double recursion 

of the course-of-values type {subscript k2). q x is the row index and q^ 

is the column index. Actually, there is no recursion on q,. However, 
since all values of b. are to be stored, advantage is taken of this property 
of the course-of-valuls double recursion. The storage index, r , defines 
the structure of the b-Cq.jqu) array as rectangular. The initial values, 

D (q ,0) — 1, are defined by means of a branch equation. The branch condition 

*f a -1* means *a < -i 1,» where a, is an integer-valued variable identified 
11 5 ~ J 7 
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with qu in the Computer Table. Since h g on the left of the equation is 
really ^(^4- 1, qg+ 1)* the condition q^ -1 is equivalent to 0^4- 1^0, 
i.e. b p (q,4> 1,0). This also explains why the lower bound of q, is -1, and 
why the index of a, is r_ ss q^ •+• 1. 

The master phase equation causes b to be computed first. It next 
refers to b.,, which leads to another phase equation. The b, -equation is 
introduced to permit the quantifier of q. to be applied to^. This is 
necessary since b. is defined by a simple recursion equation, and Rule 11.1 
prohibits more than one quantifier in a simple recursion equation. The 
output symbol, d_,_, will cause the successive values of VvSu)' *° ^ punched 
out five per cara, as they are computed. 

The above formulation is economical of time but wasteful of storage, 
since all values of b„ are stored. A better formulation, given below, 
causes b and b to be computed in one double recursion simultaneously, and 
in such a way tnat only one row of values is stored for b, and b . However, 
to obtain output, the final row of b, values must be called out In a special 
b-equation, b. • 

Formulation 2 . 
Computer Table . a Q l. Sgl, 3^50. t ^^ <Lj?0. q^O., 

r-equations . t q ~ "^j 1 * > r i = * 2 - a Q > 

^=+1.^ , \ s -l-lqg, r 5 « + la 2 , 

b-equations . b Q S. fja^y^ , 

V -i.v 5 (V-i.Vo& 2 =52 Vu\ "*' f i ' f iV 

b l -32 VllV*' f l ' + V3V1VV J > 

Voq 5 r 1 d 15 b 4 s f^, 

Remarks! To permit the simultaneous computation of b. and b , the initial 
values, b_(0,qp-4- 1), are taken to be zero. This forces us to use a triple 
branch fo? b^, necessitating tha branch equation for b,. Note that this 
equation is written outside the parentheses because it*Vill be referred to 
in the course of programming bp. m the Computer Table, one sees that there 
are two integer-valued variables, namely, a, « q^ and a^ = q^. Note also 

that a maximum of fifty data points, a^, is permitted. Since a Q and a.^ are 

constants, they each have one datum. The mY l mwn number of values of q^ is 
50. The maximum number of values of qg is kO. 
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EXAMPLE 2 . 

Given the polynomial approximation, 

arctan x s p(x) s E ' __ n A„ x , 

c^— g^ 

valid for —1 £x£l, compute a table of arctan x for — 9^ xi 9, in 
steps of 0.1. 

We introduce ADES notation as follows. Let a Q — 0.1 J b^st *, 
a-g^s A . a 2 = «/2 ; b^s arctan x; b 2 = p(x); b^= x when \x\< 1, 

and 1 /x when |xj > i; h^a 1 when |x| < 1 and -1 when |x| > i; 
Ksit/2 when x > 1, when jxj < 1, and "*/2 when x < -lj clearly, 
arctan x s b + b^pCb ). Letting f^ denote ' < S * u * i '» f gS - f !> 
f the absolute value, we have 





FORMULATION 


Computer Table, 


a Q l agl aj8 ; *jltfkfti 


r— equations . 


r^-Z.-6q. v r 2 = fj., 


b-equations . 


b o = f oV 


VoqglSo. 


d 22 \ = + V V 2 ' 




b 5 =VlO b 6- 1 " f 6V f lV 




^VllV" f i°" f iV 




^6= -"A 9 " 




b 8 - f 9 b 6' 




\ = VllV* V' f 6 x ' 




b 3 = VllV' f l b 6' A - b 6' 




b, S «bJb,, 




10 3 3' 




b ~™ ■ • b—b_ • 




D 2 * 3 9' 



Vo^eb « u + • VAoVi' f i*i r 2' 
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EXAMPLE 5. 

The usual •differential corrections » problem can "be stated as 
follows, g is a function of the independent variable, a,, and the 
m +1 parameters, ^ (0), ..., b^ (m). A table of values, a 2 (q x ) VS. 

a-ifa-j)* £ q, £ n, is given; i.e. 

a^) s -g (ajfa^), b 1Q (0), ..., b 10 (m)) . 

It is required to determine the •best 1 values of the parameters. 

The problem is linearized by writing 

a 2 (q 1 ) - B<a 1 (q 1 ), b 1() (0,0),..., b 1Q (m,0))s: s£ _ ^'AD^fi), 

where g. denotes the partial of g with respect to b^(i). The •corrections 1 , 

A b in (i), are then determined by the method of least squares, and the 

process is iterated until the difference in successive sums of squares of 
residuals is sufficiently small. One possible order of computation is as 
follows . 

1.) Compute values of b 1Q (i), starting with initial guesses, a*(i). 

2.) Compute the augmented * observational* matrix, 
< 3£ m, b^i,,}) = 6^(^(1), b 10 (0),..„ b 1Q (m)) , 

b^m* 1) = a 2 (i) - gfa^i), b 10 (0),..., \ Q fa)) , 

for O^iin. 

3.) Compute the augmented •normal » matrix, bg, by multiplying the 
matrix h. by its transpose, at the same time computing the augmented 
triangular, •Choleski* matrix, b,. 

k,) Solve the linear system of normal equations for bg SsAb.^ by 
the Cholesfci method. 

5.) Test whether the difference in successive sums of the squares of 
the residuals is less than some constant, a^. Repeat 1 - 5 if necessary. 

To write an ADES formulation, we denote the ^ and a 2 - g by the 
indexed function, f^r^, oir^m+lj m by a^; n by a ? . 
Let g s e^ + e 5 ^, where A and B are denoted by b 1() (0), b 1() (l), respect- 
ively. Thus, g s 6g /^A sxe** and g^* x 2 e Bxf ^, We denote the exponential 
function by f Q . 
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FORMULATION 
Computer Table , a^l a 1 a^l a.^pto a^O a.^', ^^a^a^ag^a^ 

r 12 9 V 50 r 9 55 r 13 9 ' 

Auxiliary Equations , (f ^ 35 QQ • c^ • c^Cg, 

f 50=09* C l' C l f e*' c l c l c 3> 

f 50 * 09 "** +f e *' C l C 2 f e ' C 5 * C lV ) 
f 51 s+Cg+l^+lc^, 

f 52 * +c 2 +1 •+.c 1 l* e 1 2* 

r-«equations . r 2 s — la^, r a, 4-lqg, 

r 5 5 + lq^ s + 1^, r ? =: f^O, r 8 a ^1, 

r 10* * lq 5' r U s + lq 6' r l^ = + ^M 
r 15 = + ley r l6 * -a^r^, r^ * -a^, 

r 12 * f lV r 15 * f lV P * B V? Vl*' 
r 9 - V^V- r l = ^^l 10 ^' 
r 19 = +1 V r l8= f lV 
b-equations . V-lVl{^" lq 2 r 2 (^lO^ = 05 P * f loV' 

f l a 3 r 3' +b 10 r 5 Vy> 
V -1^ 5 ( V-lq^a^b^ = 03 V 5 a l r 6 b 10 r 7 b 10 r 8 a 2 r 6 ,j; 

b 3 r 9 =03 P ^ f 10 a 9 a 8' /HVu. r ll? V V J J 

b U r l8 ' : - Wl/lV } 
^^VloV' V* -*l^ r i9 b i^i' 

■. \= -*2 r 10 r llV 

y**tf6 b 5 ^ll+^fu^VlO^^' *1 ' 
Vog-xo^ *ii +t 9 < iio*Vi7 r i6 b 8^io' f i°> 

Voa 7 a 5 %= 11 +^3%-^ioV > l r llV f 3°' 

EXAMPLE U. 

Given the system of n ordinary first order differential equations, 

"os-n < m, -|£- s^ (x,y Q , y r .-, y m ), 
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and the boundary values, y.(x ), it is required to obtain a solution 
by the Runge-Kutta fourth-order method. 

Using the mesh size, h, the recursion formula can be written as 
follows . 

* n£ m Ayjj + 1) * h(k ln + 21^+ 2*^ + \ n )A, 

y n (,J+ l)sy n (j)+Ay n , 

oi i£ m ^0 + 1) s'gjtoJ), y (J), ..., y m O)), 

k 2i (j4 1) S g^H- h/2, y o (j)+ 1*^/2, ...,y n (j)+hk lm /2) 
kjiU+l) CgiCxj + h/a, 700)4 1*^/2,..., yJjJ + hk^/2), 

\ ± u+ x > * %C x j + h » y ^)+ Kk 5 0'—'yn^)"*" **>)» 

To illustrate an ADES formulation, it is sufficient to consider the 
case is 2, where g Q (x,y o ,y 1 ) 35 y^x/eixx y Q , and gjfoy^y^ a y^* x q . 

The computation of y and y, will be formulated as a triple recursion in 

which the main recursion index is J . For each value of J, we compute Ay > / w « 
for all n. Then, by a double recursion, we compute the k*s. n 

The following ADES notation is adopted I 

a Q = m ; a 1 =. h. a^ ss j^(* 6 ).; ', a^s xj 0^= 3°, q^=. n ; 

b L = A y; b 2 * y,*' \q^% ~ k q5i , *^ a ^ (0) = °' *^ X) ~ 1/2 ' 
a^(2) a 1/2, & k (3) - 1,' f^ s 8i ^ = Xj . 

FORMULATION 
Computer Table . a 6 la Q l a^l a_l a^k a g 3 . a g^ a^q, r^ r_^5 r 9, 

Auxiliary Equations , (f = Q _ / - c c g sin c 1 , 

f 5o-G9 + ' c i c 2 c y ) 

f 5L s+c 8' e l + °3 1 ' 
r-^quations . ^ s f^, r ? s f^q^, r Q = 4-lq^ r Q = f^, 

b-equations . b = ^ ^n * 

V-lq 1 a 6 { VOq Q a o (b 1 r 1 s^f^O, fjO, 

• a,/ + b^lqg + • 2b 4 2qg + • 2b^3^b^6, 

V2 r i =o 5 Vio a 5 ' f i a 2 V +WlV) 
d 2ib =o 5 - , - a 3' 4 - a 5 1 V 
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Vo^C^ =03 +t) 2 q 5" W^f-k.%'} 
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